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Abstract 


A theory is formulated for the finite deformation of a thin membrane composed 
of homogeneous elastic material which is isotropic in its undeformed state. The 
theory is then extended to the case of a small deformation superposed on a known 
finite deformation of the membrane. As an example, small deformations of a circular 
cylindrical tube which has been subjected to a finite homogeneous extension and 
inflation are considered and the equations governing these small deformations are 
obtained for an incompressible material. By means of a static analysis the stability 
of cylindrically symmetric modes for the inflated and extended cylinder with fixed 
ends is determined and the results are verified by a dynamic analysis. The stability 
is considered in detail for a Mooney material. Methods are developed to obtain the 
natural frequencies for axially symmetric free vibrations of the extended and inflated 
cylindrical membrane. Some of the lower natural frequencies are calculated for a 
Mooney material and the methods are compared. 


1. Introduction 
The theory of the finite deformation of an elastic solid has been considered 
by many authors (see for example [J], [2] and [3]), and is reviewed briefly in 
Section 2 of this paper. The application of this theory to particular problems is 
difficult because of the non-linearity of the equations involved. This difficulty 
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has been circumvented in problems which have simple geometries by the use of 
semi-inverse methods, see [/] for example. = 

If the problem is such that one or more of the spatial variables can be elimi- 
nated, much simplification results. In problems involving thin shells, one of the 
spatial variables is removed if an approximate theory is developed in terms of 
the resultant forces and moments across a section of the shell and in terms of 
the deformation of the middle surface. This approach is well established for 
problems involving small strains (see [1] or [4]). A general approach to the 
exact theory of the finite deformation of shells has been given recently [4]. 

A further simplification is achieved if the shell is considered as a membrane. 
The membrane approximation assumes that the forces or stresses in the middle 
surface are large compared with the transverse stresses. The resulting theory 
involves stress resultants acting on sections normal to the middle surface and 
neglects bending moments and shear forces, and in the classical treatment of 
shells is often used to provide a particular integral from which a more exact 
solution can be derived. The finite deformation of a thin membrane has been 
discussed in [2] and several problems of this type have been solved [6]. 

In Section 3 of this paper, the general theory of the finite deformation of a 
thin membrane composed of an isotropic elastic material is re-formulated. The 
resulting equations are non-linear and are readily solvable only for problems 
with a high degree of symmetry. One means of overcoming the difficulties, which 
has application to certain types of finite deformations and also to the investi- 
gation of the stability and natural frequencies associated with a given finite 
deformation, is the theory of a small deformation superposed on a known finite 
deformation. The general theory of a small deformation superposed on a known 
finite deformation has been developed in [7]. In Section 4 the corresponding 
theory for thin membranes is developed. 

In Section 5, the theory of Sections 3 and 4 is applied to a circular cylindrical 
membrane which is subjected to a finite homogeneous inflation and axial extension 
characterized by extension ratios A, and A,. The equations of motion for a small 
additional deformation are obtained for an incompressible material. 

In the remainder of the paper, the stability and natural frequencies of the 
free vibrations of the circular cylindrical membrane with the finite deformation 
Ay, 4, are considered. The validity of static methods in stability problems has 
been discussed by ZIEGLER [8] for linear mechanical systems with a finite number 
of degrees of freedom and for linear elastic systems. In many cases a static 
method is much easier to apply than the usual dynamic method and can be used 
to determine the stability of the system [8]. In Section 6 we discuss a static 
method for the determination of the stability of the cylindrical tube in the state 
of deformation A,, Ay. 

According to the classical theory, the tube is unstable under external pressure. 
For a small internal pressure the lowest free vibration frequencies are those 
associated with highly non-symmetric modes. However, as the internal pressure 
increases, the degree of asymmetry associated with the mode of lowest frequency 
decreases [9]. Thus an analysis of the axially symmetric stability problem and 
the axially symmetric vibrations should give correct stability results and the 
lowest natural frequencies when the internal pressure is sufficiently large. 
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In Section 7 we analyse the dependence of the stability of the tube on the 
extension ratios 4, and A, for symmetric modes using a static stability analysis. 
In Section 8 we discuss the geometry of the stable regions predicted in Section 7 
for the particular case of a tube composed of Mooney material. Finally in Section 9, 
by considering the natural frequencies of the tube, it is verified that as far as 
axially symmetric modes are concerned the stable region is that predicted by the 
static approach of Section 7. 

The actual numerical computation of the natural frequencies is laborious and 
in Section 10 an approximate method for calculating the frequencies is given. 
In Section 11 numerical results are discussed in order to appraise the adequacy 
of the approximate method of Section 10. 


2. The Basic Equations of Finite Elasticity 

We consider a homogeneous elastic body which is isotropic in its undeformed 
state. 

At time t=0, the body is undeformed and the coordinates of a typical point 
in a fixed rectangular Cartesian coordinate system % are (*,, %, %3). If (9, 0, 3°) 
are the coordinates of this point in an arbitrary curvilinear coordinate system #? 
fixed in the body, we have Re gg (2.1) 


We use the convention that Latin indices take the values 1, 2 and 3 and that 
repeated Latin indices denote summation over the values 1, 2 and 3. If r is 
the position vector from the origin of the system x to a typical point of the 
body, the covariant base vectors of the system # in the undeformed body are 


or 
Ties he (2.2) 


a comma followed by k denoting differentiation with respect to #*. The element 
of length ds is given by ies pale ¢,,40' agi 


where g,;, the covariant components of the metric tensor of the system # at /=0, 
are given by 


8:7 = Ii G;- (2.3) 
The contravariant components of the metric tensor are given by 


ij ey 
Saad g = |8i; 


, (2.4) 


where c’’ is the cofactor of g;; in the determinant |g;,|. 


At time ¢, the body has become deformed and the coordinates of a typical 
point of the body in the system » are now (Xj, X,, X3). The curvilinear system 
# moves with the body as it is deformed and we have 


X, =X, (9, 0, 0, 2). (2.5) 


In terms of R, the position vector from the origin of x to the point (#, 3?, 3%) 
at time ¢, the covariant base vectors of the system # are 


GjiRy. (2.6) 
19* 


276 A. H. CorNELIUSSEN & R. T. SHIELD: 


The element of length dS is given by 
dS? =dX,dX,=G;,d0 aH, 


where 


are the covariant components of the metric tensor of the system # at time ¢. 
Finally, the contravariant components of the metric tensor are 


eS whe? 
Gia, ; (2.8) 


G=|G 


aj 


where C’? is the cofactor of G,; in the determinant |G; ;|. 

The contravariant components of stress rt’? in the system # are related to 
the stress vector #, acting on the surface #’= constant at time ¢ and measured 
per unit area of the deformed body by 


qk 


ib; G,, (2.9) 


oo VGii 


4 not summed. 
The strain-energy function W for a homogeneous isotropic elastic solid is a 
function of the three strain invariants J,, 7, and Js, 


WW. tos) 
the invariants being given by 


I, pa ion Ge 
I,=8,,G'* Is, (2.10) 
R=aGie 


The stresses t'’ are now given in terms of W and the metric tensors by 


ti = @gi+ Di + PG", (2.11) 
where 
pees ie) 
21, rai 
al, | 
Di) 
= cave 
Pe ole al,’ 


Dt i pre pirgikG 


For an incompressible material, J,=41 and 
W= Wl, 1). (2.13) 


The stress-strain law is again (2.11) where now P is a scalar invariant function 
of position and time which corresponds to a hydrostatic pressure. For certain 
rubber-like materials a strain-energy function of the Mooney form 


W=C, (3) +, (2, — 3), (2.14) 
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where C, and C, are positive constants, gives results which check well with 
experiment [/0]. The simplest incompressible strain-energy function is that of 
a neo-Hookean material, 


W=C,( —3). (2:15) 
The equations of motion in the absence of body forces are 
v1, = ef, (2.16) 


where 9 is the density of the deformed body, /’ are the components of the accel- 
eration vector in the system #, and the vertical line in ea ,; denotes covariant 
differentiation in the system # of the deformed body. If surface forces p=? G; 
act on the surface of the deformed body and if n=n; G" is the unit normal to 
the surface, where G'=G" G,, the relation 


tin, =p! (2.47) 
must hold on the surface. 


3. The Thin Membrane 


Let (#, 0?) be the coordinates of a general curvilinear coordinate system 
defined on a surface in space. If @ is the position vector from the origin of a 
fixed rectangular coordinate system x to a point (J, 9?) on this surface, the 


base vectors on the surface are 


da 
Oe An ee Oe (3.1) 


We use the convention that Greek indices take the values 1 and 2 and that 
repeated Greek indices denote summation over the values 1 and 2. 

We define a unit vector a, (91, 0?) perpendicular to the surface at each point 
in the direction of a, xa,. The distance from the point (i, 3?) on the surface 
to a point on the normal to the surface at that point is denoted by #%, and is 
considered positive in the direction of a, and negative in the direction opposite 
to a,. Then the position vectors r defined by 


r = a(9, 92) + 7 ay (9, 82), (3.2) 
— £h(H, 0) < B= Fh (H, 9), 


define the points of a shell of thickness # with #°=0 defining the middle surface 
of the shell. 

The coordinates of the convected coordinate system # of Section 2 are taken 
ho be: (75, 07, 7). 

Using equations (2.3) and (2.4), and neglecting terms in #* for a thin shell, 
we obtain 


where 


Bup—4ap, Sas—0, 833—1, 
§—4, (3.3) 
pL rl ge —0, Peds 


where a,, and a*® are the components of the surface metric tensor, 
Ag p = Ay? Ag, 42=|ayal, (3.4) 


and a*® is the cofactor of a,, divided by a. 
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At time ¢ the middle surface has a new position defined by the position vectors 
A(#', 92). The covariant base vectors of the deformed middle surface are A,, 


4. & (3.5) 


At each point of the middle surface, a unit vector A, is defined, perpendicular 
to the surface at that point and having the direction of A, A,. It should be 
noted that the unit vector Ag is not in general the usual covariant base vector 
associated with the #* coordinate curve (see also the footnote on p. 279). 

If the principal radii of curvature of the deformed middle surface are large 
in comparison with the shell thickness / and if / varies sufficiently slowly with 
# and 9, it is reasonable to make the membrane assumption that the shearing 
stress in the direction of the normal to the middle surface in planes perpendicular 
to the middle surface is small enough to be neglected in comparison with the 
remaining stresses. This means that the direction defined by Ag is a principal 
direction of stress. From isotropy, this direction is also a principal direction of 
strain. Consider linear elements of material lying along the principal directions 
of strain at a point (91, 9?) on the middle surface. Two of these elements must 
lie in the tangent plane to the deformed middle surface at (, #?). Since elements 
along principal directions in the deformed state were originally perpendicular 
in the undeformed state, the element along A; must consist of points which were 
originally along a, (#1, 92). From this and the assumption that the shell is thin, 
the position vector to a point in the deformed membrane is given by 


R= A(#4, 8?) + A, (F, G?) 8 A, (G1, 8), (3.6) 
where A, is positive and is the extension ratio in the direction of the normal 
to the middle surface. The various components of the metric tensor are now 

Gyp= Ags, Gus =0, Gsg= 43, 
G=/;4; (3.7) 
ora Getto, Guat 


ee i 
where A,, and A%? are the components of the surface metric tensor at time ¢, 
Ag = A,: Ag, Ae | Aen (3.8) 
and where A*%? is the cofactor of Ax, divided by A. 


Let the physical stress resultant m, be the force acting on the surface #* = 
constant in the deformed membrane, per unit length of the middle surface at 
a given point of the middle surface. For a thin shell, 


Tinea be, 
and defining the contravariant components of the physical stress resultant n° by 


q 1 ap 
nr, = 
Oo Ve nN Ag, 


% not summed, we have 


n= — 2h oF, (3.9) 
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We then obtai 
e then obtain n°? — A,h[Dat® + WD 4 P Arh), (3.10) 


where 
D*? = Ag a*? + {a%? at? — at" aP 1 A,,,. 


In terms of surface quantities and A, the strain invariants are 
Tee: am? Aug + And 


I, = 22 Ala, (3.11) 
1 
73 | 


Ty = Ig|aq,A%? + 


We now make the further assumption that the normal components of the 
tractions applied to the surfaces of the membrane are small compared to the 
stresses acting on the surfaces #*=constant. This means that we may take 


733 = 0, (3.12) 


Equation (3.12) then gives a relation for 4, in terms of surface quantities, namely 


PADUA Po, (3.13) 


where 
SO. Al 
D a A, 


Using (3.13), (3.11) and (3.10), we can express *? in terms of surface quan- 
tities and derivatives of the strain energy function. 
For an incompressible material J; 1 implies 
Az = (a/A)}. (3.14) 
We thus find that 
eT PAs +alA, 


| 
Vo dzp APE Ala. Ot) 


Equation (3.12) then determines the pressure function P. The stress resultants 
for an incompressible material are now given by 
n*? — (a/A8)th[{A a%® — a A**} @ + {AD*® —a A*P DIY], (3.16) 
where ‘ 
D*? — oe ar? +. f9%P gh? — gt gh 4 


fey? 


Ss pe eae 

Consider an element of the deformed membrane at time ¢ bounded by the 
surfaces #!—constant, 31+d#1=constant, #9?=constant, #?+d*%—constant. 
Let p be the resultant force acting on the lateral surfaces of the element per 


unit area of the deformed middle surface. We can write 


p=p'A,. (3.17)* 
If we also write Ae f A,, (3.18)* 


* Note that p/ and /* are not the contravariant components of p and A in the 
system #. 


280 A. H. CorRNELIUSSEN & R. T. SHIELD: 


where a dot denotes differentiation with respect to time, the equations of motion 
for a thin shell become [7], 


n°?|, + p? =h(alA)} oof”, 
n*? By, + p> =h(alA)! oof’, 


(3.19) 


where 9, is the mass per unit volume of the undeformed membrane. In equations 


(3.19), 


Bug =A * Ag p- (Ay x Ae) (3.20) 

and 
nF = OE + Dl + Tint (3.24) 

with 


Tf, = 4 A" [Ajpy + Aay,p — Apy,al- 


In a given problem, we would usually choose the curvilinear coordinate system 
on the middle surface with reference to the initial geometry of the middle surface. 
The components a,,, a%°, and a are then known quantities. The rectangular 
components X, (9, 9?) of the position vector A or a geometrically suitable com- 
bination of the X,’s may be taken as unknowns. If we write 


A, = Aie;, (3.22) 


where e; is the unit vector in the x, direction of the rectangular coordinate 
system x, and note that 
1 
A,= a [A, xy}, 
the quantities A; are readily expressible in terms of X, (9, 82). We can then 


express the acceleration components /* in terms of X, and their derivatives by 
solving the equations 


Kem fA (3.23) 


Equations (3.19) are thus three equations for the determination of the three 


coordinates X,, in conjunction with the appropriate conditions at the boundary 
of the membrane. 


4. Small Deformation Superposed on Known Finite Deformation 
of a Thin Membrane 


For a thin membrane, we consider a state of deformation which at each 
instant of time differs only slightly from that of a known finite deformation. 
At time ¢#, the vector from the origin of the rectangular coordinate system %* 
to a point on the deformed middle surface is written as A-+¢A’ where A is 
the position vector of the deformed middle surface for the known finite defor- 
mation and ¢ is a constant which is considered small in the sense that second 
and higher powers of ¢ may be neglected in comparison with e. It is convenient 
to represent A’ as 


A’ = 0, A’, (4.1) 
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where A*?= A, (w, are not the covariant components of A’). The covariant com- 
ponents of the deformed surface metric tensor are now A,,+¢Aj, and in terms 
of w,, 


Aug = Ws,p + Wpa— 2 Kb pa Bh (4.2) 


The contravariant components of the surface tensor and the determinant of 
AypteAyep are now A*®+4¢4'*8 and A+eA’ respectively. We find that 


A’11— A1[A,.—A’A,,/A], 
A’12— A4(A‘A,,/A — Aye}, (4.3) 
A’22 — A4[Aj,— A’ A,,/A], 
and 
A Va Aes + Aig tee 2 Ais: 
The strain invariants are I, + ¢I;,, where 
Ty = a" Aug + 2Ag As, 
Ty = (A,/a)(24gA +A, A'], (4.4) 
T= I,{a,p4'*?— 28} + e{agg A +35}, 
the extension ratio in the direction perpendicular to the middle surface being 


Aste dg. Writing the stresses as t'/+ et'", 


pea OD a yep Cpt ny Dy eh pA’ ee (4.5) 
where 


r , , / PD / 
QD’ = 215° (WintWels aie W, 31s] Bear aht 


Lt / / / vy Le 
Ne Ae eg EE era A iS ceeioyahe) (4.6) 


1 , / / 1 / 
P’ = 214 [Wg li + Weala + Was ls] +57-Js. 
3 


with the notation 
ow 
i de . 
v7} él; al; (hh, 2» Is) 
Also, 


D’ apes OAnAs xh dis {76 gt” — axl af” Ae. (4.7) 


The equation 133+ ¢7’8=0 gives the following relation for As, 


, As 
@’ ce Co Ag +3, P+ av? A’, — of O. (4.8) 
3 


a 


. / 
The contravariant components of the stress resultants are nF en’ *® where 


ni’? — ji(ajAyila’e? 4 4 gre], (4.9) 
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For an incompressible material, 


, 1 a \% 
Pedals sill fassicie sy Nik 10 
i= — 2 (aa) 4 (4.10) 
In this case, without knowing the explicit form of the strain-energy function, 
we can eliminate 43 and P’ to obtain 


veh oat 6 ane] 4p D*? — AGP ab? AL.) + 

(4.11) 
a v4 ys [4 sak — grab tap @ yap mv g! 

+4. (6+ a Ay, FI 4 A |+¥|D Vr ace 


where 


@' =O, + @,, 12, 


; / / oP 
ve = ©,,1,+ PJs, Oy = ae Ua 1s) - 


(4.12) 


Ti, and I are now defined by (4.4) with A, and A; given by (3.14) and (4.10). 

The force vector acting on the lateral surfaces of the membrane is now 
(p+ ep'") (4,+ Aj) and the acceleration vector is (/*+ef'*) (A,+ Aj), where 
Aj,=A',. We find that the equations of motion for the small additional defor- 
mation are 


‘ap , , , , , 
On ats fg AO De go te gh Eg ea 


Col ad 
2 , 1 A’ 
= hol) |f°->4F,. (4.13) 
a , lo , 3 , Ay 
n Bip tn t Bee + p *= he ja s— 5 =P, 
where 
De = A: A ies ate Aas re Aby,al at LN oe + Ajy,6 x Ag,,a] , (4 14) 
, 1 / , ’ A’ : 
Ep [A.p° (Ay X Ay) + Ay pg: (Ay X Ae +Ay x A,)]— Sa Pas: 
Equation (3.23) becomes 
Kot SS ee ee 
and at once we have 
Xi=f A+ Pf Ag’ (4.15) 


a relation to determine /’*. We can easily obtain Aj‘ using the relation 
/ | , ’ A’ 


Expressing X; in terms of w,, we obtain from (4.13) three equations for the 
determination of the three components w,. 


5. The Circular Cylindrical Membrane 
We now consider a thin membrane which in the undeformed state has the 
shape of a circular cylindrical tube of length J of constant thickness. Let the 
radius of the middle surface be c and the thickness of the tube be h. We will 
choose the right-handed rectangular Cartesian coordinate system x so that the 
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axes 1 and 2 lie in the plane of one of the ends of the tube while the third axis 
lies along the axis of the tube in the direction of the other end of the tube, 
Fig. 1. The convected coordinates of the middle surface are chosen in the follow- 
ing way. Let # be arc length measured from the plane x,=0 along the line 
formed by the intersection of the middle surface and the plane x,;=constant. 
0 is the distance from the plane x,;=0 measured in the direction of the x axis. 

The position vector a to the 
undeformed middle surface is then 


; as 
a= ccos—e,+ 
Cc 


(5.1) 


at 
+ esin Z e, + #e;, 
C 
from which it follows that 


1,=1, A&o=0, 4.=—1, 


Ae a2—Q, q22— 4, (5.2) 


x2 


a=1. Fig. 1. Coordinate systems for the circular cylindrical membrane 


We consider the quasistatic homogeneous finite deformation of the tube by 
a simultaneous inflation and extension. The radius of the middle surface becomes 
A,c, the length A,/, and a point originally with axial coordinate x; now has 
coordinate X,;=A,x3;, 2, and A, being positive constants. We may write the 
position vector to a point on the deformed middle surface as 
ot F ot 
A = hc cos — @, + Ay c Sin —— €y + A, OP es. (5.3) 
It follows directly that 
Ay=4M, A;,=0, A,,= 43, 
Al =, Alt— 0 pA (5.4) 
1 
AA 2. 


For the strain invariants J,, we find that 
p= A+ A+ 43, 
Ty = Aj Az + AAS + ALAS, (5.5) 
Ty = Aj AR AG, 
where A, is given in terms of A, and A, by 
B+ (i+ B)Y +5, P=0. (5.6) 
The stress resultants are then 
mt = Ah | D+ (+ 23) ¥ +55 PI, 
WOE 


nt = Mh |®+ (A + 48) Y + zy P|. 
iz 
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Also we have 


Ba = ’ Big Beg 0- 


The forces necessary to support this deformation are a normal pressure on 
the lateral surfaces of the tube given by 


path) (g++ v+s, Pl, (5.8) 
and an axial force T per unit length on the ends of the tube, 
T = 1h yh] ® + U8-+ 28) ¥ +35 PI. (5.9) 
2 


For an incompressible material, 


L=A+4+ ; 


nial (5.10) 
1 1 
The stress resultants become 
m= 5h [a 1)[6+ BY], 
1 
n2—(, (5.41) 
nes — Wai aa ater aT. 


The forces necessary to support the deformation for incompressible materials 
are given by 


p= ay MA 118+ BY), (5.12) 
io aE [az as — 4] [@ + 22]. (5.13) 


We will now give the steps leading to the equations of motion for a general 
small deformation superposed on the above finite deformation. It is convenient 
to express the components w,, defined earlier by A’=w, A”, in terms of new 
dependent variables u,v, and w. If a cylindrical polar coordinate system is 
defined in the usual way with respect to the system x, ew, ev, and ew are the 
physical components of displacement (in the polar system) of the small additional 
deformation, with w radial, v azimuthal, and w longitudinal. We find that 


2, =Ag?s 


W,=A,w, (5.14) 


Ws =U. 
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The perturbed components of the deformed metric tensor are from equations 
(4.2), (4.3), and (5.14), 
Ain=2h(S+0i),  Ata=Ay.+dewr, 


Az, = 2A, », 
2{[u | (5.15) 
cA oes eee Es fly ee 
ks a Uy A 23 22 [Av otApwil, 
, 2, 
A oe ane 2° 


and 
AA A la. “ + Ayo, + A, w,2], 


where v , denotes dv/08", etc. For Ij,, we find that 


T= 244s + 24 [% +0, +420, 


Ig= 2h Whey Ma + Ao{% +014 wd, 
(5.16) 


Ig = 24 (A§ + 28) w+ 22,03 +28) w, 
} 28) 0,1 + 225 a (AE + 28), 


+ 2A, (A434 


where we may evaluate Aj in terms of u,v, w, and known quantities by using 
(4.6), (4.8), (5.2), (5.4), (5.6), (5.15), and (5.16). We also find that 


DU=A4+43, Dit=0, DB=AiTA;, 


D'1 = 2/313 + 22,09, 
(5.17) 
DAs RN Ay 1, 
, , U 
D 22 — 25 Ag+ 2h (= +2). 
From (4.5) and (4.9), we can now calculate the components n'*8. From 


equations (4.14), 


ul 1 Uy ip (Poe 1 Uu ] 
= | ee Vv 5 => — | —— Uv 
11 + | pe Be gah Sapa ae 12| 
al 1 es 1 Ay “4 
9 = UV ) Ww 
Ids 7, 2? 11 | ge eee 
5.18) 
1 2 (5 
Tis a ia: 132 = | #22: 
2) 
Bu=— C11 — 20, — 2, Bye =— [CU y2—% ol, 
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Finally, from the definition of the finite deformation and from equation (4.15), 


we obtain 
pao, ft=td, ptata,  fi=i. (5.19) 
2 


Substituting into equations (4.13), we obtain the equations of motion for w, v, 


and w. 
For an incompressible material, we find that 


, 1 | 
A3 = — 22 [aa(% +0,] + A, w 9] « (5.20) 
The quantities J’, can now be written 
' Ai AZ—1 [a Ai Ag—1 
h=2 Ae fe tal +2 gg ee 
741 274-1 oe 
, == Uu Ay Ag 
Nine PD 2B E v4 2 z 72 Ws 
The components n’** become 
Y, h y. 7 | 
nit = ee [5 — AA) (® + BYP) + 2K,][~ +0,| + 
142 € 
h 
Su We [(3 — At Ab) ® + (ALAS + 1) AY + 2] Wo, 
Ne ee an [D+ (Ai + As — At 22) )¥] [Ay si pet Aw il, (5.22) 
/ h 2 a 
m8 = sare [5 — At Aa) (P+ AY) + 2Ky] w, 9 + 
i 2 ‘ 3 
+ ayy (OHA) © + (Has +1) RY +2H)[* +0], 
where 
[Az As —1] [At AR— ‘ 9 Sane 
HA 1, + (+) ®, + BBY], 
1 22 22 [ I, sine 1 A Ty | > (5.23) 
442 475 
K,= [Ay AZ—1]? [®,. + 243 @,, + 18%]. 


22 72 


After some reduction, the equations for small vibrations in the case of an in- 
compressible material are found to be 


By W 22+ : Ust B12 + Pg Mir +— - a Pp? =A A ew, 


Bs 99 — ie Ce By W 9 — fs v i +Beu, Lact nue pie a 20%, (5.24) 


B; Wi2 + Bs ( + 011)+ Bs, 22+ ae pra he 205, 


where 
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By = (3 + ALA) (G+ AV) + 2K, 

Bo= 7? (3 — AtA2) ® + (043-441) BY + 2H], 
1 

Bs= (AAS — 1) (B+ HY), 


2 
Be YP KG, 
‘i 
re (5.25) 
Bs= 59 (B+ Ai Aa) (P + AP) + 2K), 
Ba = 72 (a8 8 — 1) (G+ BY), 
Mal 
p= 38+ (+ + BAY + 2H), 


Bs = Aa (Ai AO +), 


and where @= 9, is the mass per unit volume of the material. 


The constants £, to Bg can be written in terms of the pressure #% and the 
axial force T of the finite deformation as follows: 


A RB AB c ap 
Mitiety s Cie 


3°43 
Be wv iT 
gee de® 
: oR ' 
B= Aj AR c (A, p*) 
h GIN 


AARC 


al 
WO 
= 


Be= h P, 
Feieating 2 Bac ap? 
Pasehh aaesa Rinse Des ols 'asieih. eed hee 
3°95 
b.= se Gate Anak). 


h (A2— 22) 


From equations (5.26), we see that equations (5.24) could have been derived 
directly by considering the motion of an element of the shell. 

Note that f,, 8,, and B, are always positive for a Mooney material, defined 
by equations (2.14). We will assume in the following that they are positive for 
the incompressible material considered, 


By, Ba, Bg > 0. (5:27) 
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6. Stability 


Consider a body in a known (static) state of finite deformation. If we apply 
external forces to produce an additional small deformation (subject to certain 
boundary constraints), the body will in general vibrate about the position of 
the known finite deformation upon removal of the external forces. If the result- 
ing motion of the body remains bounded for all time and for every choice of the 
disturbing forces, the initial state of finite deformation is said to be stable for 
the particular constraints. 

Writing the displacement vector of the vibration at time ¢ as eR’=ew,G’, 
we may obtain the three free vibration equations for the components w,; by using 
the theory of a small deformation superposed on a known finite deformation [7]. 
The boundary conditions on these equations are determined by the constraints 
and the initial conditions are arbitrary within the class of initial conditions 
possible for these constraints. 

Writing w, = W, (9, 02, 03) givt 
and noting that the equations for the components w, contain only second order 
time derivatives, we find that w? may assume only a certain sequence of values 
win)» n=‘, 2,3,..., which depend on the particular finite deformation and on 
the form of the constraints. If all the values w/,) are positive, the finite defor- 
mation is stable for the given constraints but if any of the values w7,) is zero, 
negative, or complex, the deformation is unstable. 

A body is said to be statically stable for a given finite deformation and certain 
constraints if no non-trivial equilibrium configurations exist neighboring to the 
finite deformation and satisfying the constraints; it is statically unstable if such 
configurations do occur. 

Any state of a body which is statically unstable is unstable since there will 
then exist an integer m such that Win) =. On the other hand, a statically stable 
state may well be unstable. 

The concept of static instability is valuable since it frequently enables us to 
discuss the stability of a system without investigating the reality and sign of 
the values w/,). ZIEGLER [8] has shown that for a linear, non-dissipative, non- 
gyroscopic system with a finite number of degrees of freedom, the stability can 
be determined by a study of the non-trivial equilibrium configurations (latent 
instabilities) of a system. An example is the Euler column with one end fixed. 
A static stability analysis gives the well-known result that latent instabilities 
occur for a sequence of axial loads L;}<1,<L,.... An analysis of the vibrations 
of the column with an axial load shows that the column is stable if the axial 
load L is less than L, and unstable for all larger loads. 

We now consider the stability of the circular cylindrical membrane of Section 5, 
composed of incompressible material, for the finite deformation consisting of a 
homogeneous inflation and extension defined by 4, and 4, and for the constraint 
of fixed ends. Setting 

“= UH, 0%)? 


v= V(H, 3) ery, (6.1) 
w= Wh, 8?) eek 
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the equations for small vibrations are from equations (5.24), 


PWoot Us + ByVistPpWu=—eW, 


BoUes— U— Bw, "V+ B.Un=—MU, (6.2) 
B, W. 12+ Bo(—2+Vin1) + BeVio2=— V, 
where 
v= 0 At Aw. (6.3) 
The boundary conditions are 
Opa VO at, Oa. (6.4) 


Solving (6.2) with the boundary conditions (6.4), we obtain a sp dltaes: of values 
Vin) that v2 may assume, and corresponding eigenfunctions U™, V™, and W™. 
Defining Q as 

Q= (on) ae Vin)) (Uu™ UM a vm y™ 2h w) w) ‘ 


we find from equations (6.2) that 


ae ae (m) ow™ (n) ow™ | rs) f ion) eum (n) oU™ ' 
Cae arti Op YP age f + Ps gga \U 69 a0? \4 


0 (m) eV™ jagcy 2 0 (m) 0V™ 7 (n) CVE Ss 
Peart! ger Var} Ba ag Ga — Sart 


Bo 8 grzim) yyy(m) _ gylm) rn) eye ny BU") 
soa ogg algal pais aru ELST eae 


é (my) OW™ (n) EW) Bs 8 grrtom) yim) _ gym) gym) 
Ps ae om ag } i. Hel {J U WAC }+ 
Spee pee ey) shit { (n) W™ in) OV™ | 
Bb, Ov {v Bry T W agl Bb, ap V ape W a [- 


Integrating Q over the (deformed) middle surface of the membrane, using (6.4) 
and the fact that the eigenfunctions Uu™, vy Ww are periodic in #1, we have 


1 2ne 


f fOdPd#=0. 
0 0 
sthat-1s; 


1 2ne 


= (v7, in) Ven ) J J [u' m) [7 (n) a vm yo ue w™ Ww”) dm de. 


Since the coefficients of equations (6.2) are real, we know that if 17,,) eyo 
WwW) is a solution to equations (6.2) under the conditions (6.4), ee Ven) 5 ID: 
v™, Ww) is also a solution, a bar denoting the complex conjugate. Hence 


1 22c 


OF (%m e Yn)) ut J [| Ue |? a \v™|2 me |w™) |?] de dg. 


Since the integrand is non-negative, it follows that 17.) =n). Therefore, for the 

finite deformation and constraints under consideration, »7,,. (and therefore «/,,)) 

is real for all m. The deformation characterized by 4,, A, will then be stable if 

all the values »7,,, are positive and unstable if one or more are zero or negative. 
Arch. Rational Mech. Anal., Vol. 7 20 
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Each finite deformation characterized by the extension ratios A,, A, can be 
represented by a point ina A,-A, plane. The undeformed state is then represented 
by the point 4, =2,=1, and possible finite deformations by points in the quadrant 
A,, Ag> 0. In this plane, the stable states of deformation for fixed-end constraints 
will consist of a region of the plane including the point 2,=/,=1. 

Suppose the lines of static instability are drawn in the A,-A, plane. Consider 
the region consisting of all points which can be reached from the undeformed 
point without crossing a line of latent instabilities. This region clearly is bounded 
by lines of latent instabilities and possibly by the edges 4;=0, A,=0 of the 
quadrant. If we consider a particular value »? as a function of A, and A, the 
lines of latent instabilities are located where »? changes sign. Thus, since no 
line from the undeformed point to any point in the above described region 
crosses a latent instability, we expect that all the values 77, are positive through- 
out this region, 7.e. that the foregoing region is the region of stability. 

In Sections 7, 8, and 9 we will discuss the stability for cylindrically symmetric 
modes of deformation of the circular cylindrical membrane for the finite de- 
formation A,, A, with the constraint of fixed ends. 

In Section 7 we will obtain the stable region by considering the latent in- 
stabilities, while in Section 8 the geometry of this region will be discussed for a 
Mooney material. In Section 9 we will show that the region given by the results 
of Section 7 is indeed stable by determining the region where all values y7,) are 
positive. Further we will show that for certain values of the ratio C,/C, for a 
Mooney material (see equation (2.14)) there is another stable region separated 
from the point A, =A,=1 by lines of latent instability. 


7. Static Analysis of Stability 
In this section we determine the largest region such that every one of its 
points can be reached from the undeformed state 4; =/,=1 without crossing a 
line of static instability. 
The equations for the determination of symmetric latent instabilities are 
from (6.2), 


(7.1) 
Bs 4 o2— Bay — Bs ws =0; 
and 
BsV22=0, (7.2) 
with the boundary conditions 
“u=v’=w=0 at #=0,1. (7.3) 


Note that the equation for v does not contain terms in u and w, and vice versa. 
From our assumption (5.27), we see from (7.2) and (7.3) that the displacement v 
does not give rise to any symmetric latent instabilities, 


For equations (7.1) we assume solutions of the form 


id 
WZ BIER, WED, 
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where B,, B,, and « are constants to be determined. Substitution into equations 
(7.1) shows that the values of « for which the constants B,, B, can be non-zero 


are given by 
= 1 f Bi Bs—Be \3 
el rab cal fe 


The type of solution depends on the sign of (8, By —3)/B; Bs. 


(a) In the case when 


dn Ba Ba miPies, ola 
<0 = 0; =i 0); 
C B, Bs ” 


the solutions for wu and w are 


w= A,+ A, 0+ A;sinhy # + A,coshy 8, 


pest Bs A,—c Py y [A,sinh y 3? + A; coshy 3], 
By By 

where A,, A,, Az, and A, are arbitrary constants. A latent instability is present 

if A,, A,, As, and A, can be chosen so that the boundary conditions (7.3) are 

satisfied. The requirement that non-trivial equilibrium configurations exist is 

found to be 


g(y) =coshy —1— Bysinhy=0, (7.4) 
with 
a — BBs 
y= l eaery : 


The parameter B is positive because of the assumptions (5.27). Since 


Finf2 = 1 2B 


equation (7.4) has no real solutions for B>}. Suppose that £, 6,—f3>0 and 
Boe.0...1 hen; 


B,Bs> 62 or Baie >t. 


Therefore in the region f,> 0, 8,8, —/2> 0 there are no latent instabilities. It 
is easy to show that latent instabilities do exist in the region 


Ps <0, B, Ba — B2< 0. 
(b) In the case when 


1 ~BE—B, Bs — §2?>0 O10 
c? B, Bs 


the solutions of equations (7.1) are 
w—A,+ A, + A,sin 60% + A, cos 60, 
u= — c Fe Ag+ ot 6[A,sin 60 — A, cos 687], 
4 2 


20* 
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where A,, Ay, Ag, and A, are arbitrary constants. The condition for a non- 
trivial equilibrium configuration 1s now 


f(y) =1—cosy — Bysiny =0, (7.5) 


where 
By B 
y=0l1, B= 5 BB ; 
The initial dimensions of the membrane appear in (7.5) only in the ratio J/c. 
Thus the lines of latent instability are the same for each cylinder with the same 
initial ratio of length to radius. For a given value of //c, we see from (7.5) that 


Ol=2n7, ij Ni Gane 


defines some, but not necessarily all, of the real solutions of equation (7.5). 

Consider a point P in the A,-A, plane which lies on that part of the line 
f3=0 where f,B,—f3>0. Let Q be a point, in the region B;<0, B, By —f3> 0, 
lying arbitrarily close to P. At Q, for a given J/c, dl is positive. If m is the 
smallest integer such that 6/ at Q is less than 2mz, 61 will assume all the values 
2nx, n>m, as we move along the straight line from Q to P since as the point 
P is approached, 6l/->+ co. That is, any point on the line 6;=0 for which 
B, Bs —B3 > 0 is a limit point, for all //c, of latent instabilities in the region fB,< 0. 
However, it can be easily shown that no latent instabilities occur on the line 
f,=0 itself except where £3 —f, B,=0. 

Again let us consider a given value of the ratio J/c. We find that 


(0) = 4£ (0) =0, 


while 


a2 f : 

oe =siny[1—2B]+ Bysiny. (7.6) 
Consider the part of Region (b) where B, and £3—,f, are positive. Here, 
0<B< 3. Thus d*//dy? is positive for all y in the range 0 to z. For a< y< 2a, 
f(y)>0. Therefore, in the region where f, and £3—f, B, are both positive, the 
line of static instabilities corresponding to the smallest value of y=6/ is given by 


B3—B, Ba 2(o 
pee es 1), F 
FES = an? (cll) (7.7) 
If a small additional normal pressure ¢/’? is applied to the inner surface of 
the finitely deformed cylindrical membrane, the solution for the normal displace- 
ment w with w=0 at #?=0 and J, when the point 4,-A, is in Region (b), is 
given by 


— cos d& — cos — (7.8) 
where f(y) is defined by equation (7.5), and where 


=H 12 Bt ee ay Ok 
GS / Bh ol p sin >. 


From equation (7.8), we see an illustration of the well-known result that when- 
ever a system possesses a non-trivial equilibrium configuration GOO), the 
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application of small surface forces would cause an unbounded deformation. It 
is easy to see from (7.8) that for B;> 0 and for 


\<e BiB, By < 42? (c/l)2, 
B, Bs 
u is positive for all pertinent & (or 92). However, when 


BBP > 450 (c/l) (7.9) 
1P3 

the radial displacement ew has one or more interior nodes (where w=0), and is 
an oscillatory function of #?. This behavior of the displacement ew suggests 
buckling (and instability) of the tube wherever (7.9) is satisfied. 

It can readily be shown that for £,>0, the line B3—, 8B, = 0 has no static 
instabilities except for the limiting case //c->co where the line B3—f, B,=0 
coincides with a line of static instabilities. The intersection point of B,=0 with 

5 —B, B,=0 has a latent instability. 


Consider the region defined by 
Bs20, B2—B, Ba < 4B, Bs 2? (c/l)?. (7.10) 


From the results of this section, there are no static instabilities within this region. 
Moreover, as can be seen from equations (5.25) and (5.23), this region includes 
the point 4,=A,=1. Starting at the point A, =A,=1, we find that we can not 
leave the region defined by (7.10) without crossing a line of static instabilities. 

If region (7.10) is simply-connected, it is the region defined in the first para- 
graph of Section 7 and is the stable region, as will be shown in Section 9. If 
region (7.10) is multiply-connected, the region we set out to find is the part of 
(7.10) contiguous to 2,=/A,=1 and is stable. However, as we will see in Section 9, 
the ccmplete stable region consists of (7.10) in its entirety. 


8. Stable Regions for the Mooney Material 
In this section we will investigate the geometry of region (7.10) for the 
Mooney material, equation (2.14). The ratio of the constants entering into 
expression (2.14) for the strain-energy function will be denoted by I’, 


Ds Cle (8.1) 


The line B,;=0 is a (closed) boundary of the region (7.10) and is the same 
for all incompressible materials and for all ratios J/c. From equations (5.25), 


this curve is given by poe tO. (8.2) 


and comparison with (5.13) shows that it is the line on which the axial force T 
in the membrane is zero. It is well known that a thin membrane becomes unstable 
whenever it is stressed by compressive forces, and under axial compression the 
cylinder can buckle in symmetric modes. The cylinder is unstable also under 
external pressure, that is for negative ~3, and the cylinder then buckles in un- 
symmetric modes. From equation (5.13), the pressure is zero when 


Aha 0; (8.3) 
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independent of the form of the strain-energy function and of J/c. The present 
discussion of stability is confined to symmetrical modes in the dynamic analysis 
and to symmetric latent instabilities in the static analysis and the instability under 
external pressure does not appear. It 
must be remembered in the following, 
therefore, that points of the region 
(7.10) for which p?<0, that is below 
the line (8.3), are unstable. 

The other boundary of the stable 
region is not so obvious from physi- 
cal considerations and depends strong- 
ly on the form of the strain-energy 
function. In order to discuss the shape 
of the boundary (7.7), it is useful to 
first investigate the curve B3—f, B4=0. 
This curve depends on the particular 
form of the strain-energy function used 
but is independent of the ratio J/c. 
Fig. 2. A,—A, plane and region (7.10) for 2 Mooney material Fig. 2 shows the curve 2= O and 
with 0<I"<1. Bg=0, «+++ PF=0, Region (7.10), : 

—— B,=0, - - - - 82-8, 8,=0, —-— B2—B, Bs=4 Bi Bs 2? (C/L)? indicates the shape of the locus 

pe— Pi pa=0 for 0O-< I< 4." It tam 
be shown that the locus consists of three branches, two of them asymptotic to 
the curve 6,=0 for large A,. The curve 6;=0 has two intersections with the 
curve f3—f,B,=0. The parameter f, is greater than zero above B,=0 and 


g 

, |@ Intersection of Bf-B, 8-0 with B,=0 
(2)----- Vertical asymptore of 7-8, By =0 

7 (3)—— Intersection of BF-B, Bu= with A,=0 


=. 
“= 
T=... 


Oo 70 20 30 40 A> 
Fig. 3. I’—A, plane for Mooney material 


negative below, while £3 — B, By is positive above the upper branches of Ps —B,Bs=0 
and below the lower branch, negative elsewhere. 

B or very large values of J/c, the (open) partial boundary of region (7.10) 
defined by equation (7.7) virtually coincides with the part of the curve £5 — 8,8, =0 
oe f3=0. As l/c decreases, these boundary curves separate from the curve 
Be —B, By=0. In Fig. 2 the branches of the curve (7.7) are indicated for an inter- 
mediate value of the ratio //c. For very small J/c, the upper right-hand branch of the 
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curve (7.7) moves to the region where A, is large, the other two branches approach- 
ing the curve 6s=0. The shaded area in Fig. 2 indicates the shape of the region 
(7.10) for values of J” between zero and unity, and an intermediate value of J/c. 

In Fig. 3, a J-A, plane 
is taken and the following 
curves are plotted: 

(1) For a given J’, this 
curve gives the value of A, 
where $3—6,f,=0 inter- 2 
sects B,=0. 

(2) For a given I’, this 
curve gives the values of 2, 
where the upper branches 
of 63 — B, Ba=0 have verti- 
cal asymptotes in the /,-A, 
plane, Fig. 2. 

(3) For a given I’, this 
curve gives the value of % 7 2 3 I, Z 
4, where the lower branch gig 4. stabte region for Mooney material; '=0.1, I/e=6, 00 —— f,—0, 
of B3—f, B,=O0 intersects ~~~ p?=0, —— B2-B,Bi=4 By By 7 (6/1)? 

Bs=0. 

As I" increases from the particular value used in Fig. 2, there is a general 
shift to the left of the branches of 83—, B,=0, and the locus (7.7) still bears 
in general terms the same 
relationship to 63—f,B,=0 
as we have outlined above. 
For decreasing J’, the shift ~” 
is to the right. 

Fig. 4 shows the 4,-A, 2 
plane for the particular 
values = 0.1 [6] and d/c=6 
within a physically reason- 
able range of /, and A,. 
The intersection of the 
lower branch of the curve 
Bs — B, Bs =0 with A,=0 
occurs for a value of A, 
greaterthan8.Anextension 4% 5 2 3 a 
ratio of 8 is too large to be Fig. 5. Stable region for neo-Hookean material; //c=2, 6, co 
of physical significance. For ---- p°=0, —— B38: Bu= 4B, Bs 2° (c/l)? 
llc< 6, the upper right- 
hand branch of the curve (7.7) lies to the right of A, = 7.87 so that, unless we are 
dealing with a very thin tube, this branch is unimportant. Also shown in Fig. 4 
are the curves (7.7) for J/c=oo and the curve (8.3) for which £?=0. 

Decreasing J’ to zero, we obtain the case of a neo-Hookean material, see 
equation (2.15). In Fig. 5 the relevant curves for tubes of various initial ratios 
are shown for a neo-Hookean material. 


g 


3 


B3=0, 
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One major difference between the stable regions for the Mooney and the 
neo-Hookean materials is that for the former there is a range of values of A, 
for which the membrane is stable for all values of A, which are large enough. 

From Fig. 3, we see that for [’>8.72, the curve f63—B, B,4=0 intersects 
f,—=0 four times. This occurs because the peak of £3 —B; By =0 has moved so 
far to the left that it has crossed B,=0. These values of /’ are far larger than 
those of physical interest but it is of theoretical interest to note that in this 
range the stable region (7.10) becomes doubly-connected and we have a stable 
region separated by instabilities from the stable region about 24;=/,=1. 


9. Dynamic Analysis of Stability 
In this section we will show that the region (7.10) is the region of stability 
for cylindrically symmetric deformations by considering the sign of the frequency 
parameter »?= Aj 23 9 w?. 
The equations of motion are, from equations (6.2), 


B,W.oo+ 2 U, =— PW, 
(9.1) 


B B; 2s 
Bs U os a Clee: We =o 
and 
BeVie2=— VV. (9.2) 
The boundary conditions are again given by equations (6.4). 

The equation for the circumferential displacement V is uncoupled from those 
for the radial and axial displacements, U and W. It has been seen that this 
mode has no static instabilities. From (9.2) and (6.4), we easily find that the 
frequencies for circumferential modes are given by 


e Dee 
ene ae (03) 


and, from the assumption B,>0, we see that v2 is positive for all n. Hence the 
membrane is stable for all symmetric torsional modes. 

Let us now examine the equations for U and W. Along the line p,=0, i.e. 
for zero axial force, equations (9.1) become 


B.Woa+2U,+W=0, 


(>—£:)u- 2 w,=o. 
Assuming solutions of the form 
U=B&™ 


W = B,e*™, 0-4) 


where B,, B,, and & are constants, we find that for non-zero values of B, and 
B,, « must satisfy 
arya 7) 
ce 


i yp, 


—s (Px Ba— 63) —B,v® 
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(2 Ps) 
Cc 


[= { > 5) 
x (Bi Ba B3) —B, v? 


Setting 


the general solutions for U and W are 
U = B, sinh w 3 + By coshu 0, 
W = Bssinh w 0? + B,coshu&, 


where B, to B, are arbitrary constants. Applying the boundary conditions, we 
find that they require 


sinh ul=0O, 


and it follows that ~/=nai. The frequency parameters v2 are therefore given 
by (two values for each integer 7) 


n2 m2 n2 m2 \2 n m2 4 
(Ps 4 Pr (ic)? iz) 4 : {(Bs4 Pr (ic)? 2) | 4 ey? 0} | Oo) 
where 0 = 83 — 8, B,. For O<0, there can be instability only if for some integer m, 


2 m2 2 2 2 
(Bs + Br Giaya) <4 -Gyays (Bx Bs — 63). 
However this is equivalent to 
n2 Ch 2 n2 m2 
Be — Bs Gigs] <— 4B aae: 


which is impossible. Therefore points along the line /;=0 are stable for 


Bs —B: Bs<0. 

For B3—f,B,>0, v%, has one negative value and one positive value so the 
points on the curve £,;=0 are unstable for B3—f, B4>0. For f,=$3— 6, Bs=0, 
we have instability since from (9.5), there must be zero natural frequencies at 
this point. 

Consider a general point where f; is non-zero. Assuming solutions of the form 
(9.4) in oe (9.1), we find that « must be of the form 


Pie re [B1 Ba — B2 — (B; + Bs) 0?» + 
+ {(93 — PBs + (Bx + Ba) 9°)? — 48, By (c9* — By) cP™}!), 


For each value of v, there are four solutions for « and they will be denoted by 
Oy, , %, %,. From a dimensional argument we see that both »* and a are 
proportional to 1/c?. Defining 6, and 06, as 


2 
Vy 


oe 


(9.6) 


= 0;,79 G2—=— O14; %=9,, %=— dy, (9.10) 
we find that the general solution of equations (9.1) is 


CA A em TERIA eMC Ae Tn? 
[Oa TAS 3 (0.41) 


— 6,9) 
—Ay,e So 


2 


fe LM pe ek DN CRY i eae 


58 
G By O24 9? L (Ang 


By aoe 
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where A, to A, are arbitrary constants. In order that the fixed-end conditions 
be satisfied, we must have the following relation connecting v®, 6,, 62, and the p’s. 


0 = [63 (B, 53 + v2)? + 68 (B, 6% + »®)?] [e*? — e %7] [eo %! — eb] + 
+ 26, 55 (By dt + ”) (Bs 68 +9”) x (9.12) 
x [fen#! — oN fehl fal 4 feb! __ gil) fof! eh, 
This can be simplified to 
(6, + 65) (v? + B, 6, 2) sinh a ; es + (6, — 65) (v? — B, 6, dg) sinh = 2}. (9.13) 


Equations (9.6) and (9.12) or (9.13) completely determine »? and in (9.12) 
and (9.13), / and c occur only in the combination J/c. . 

First we will show that the region f,;<0 is unstable. Assume that in this 
region there is a solution vy? which is negative. Then, as can be seen from (9.6), 
we can write i= bn eae, 
where 6 and d are real and positive. Equation (9.12) becomes 

Sia) oa 6, 
S (v?) = [d? (v2 + B, 67)? — 62 (v? — B, d?)?] sinh 6/ sin dl + 
+ 2d d(v? + B, 6?) (v? — B, d?) [cosh 6] cos dl — 1]. 


where 


For |v?| very large, we find from (9.6) that 6? and d? are given by 


2 AE py 2 
y By P By (Bi:— Bs) c2 (By —B3)* v2 c4 [B2 — (8, — Bs) Ba] + 


+ gaps 1B — (Bs — Bs) Be] [B8(Bs + Bs) — Bi Bu(Bs — A). 


lee : 2 BS 2 
d By BBB) [(B: — Bs) Bs — Be] (B—B,)2 vc [B2 — (B, — Bs) Ba] + 


3 


apiece 88 — (Pr — Bs) Bal (82 (Bs + Bs) — Bi Ba lB — Bad 


where terms involving »~® and higher negative powers are neglected. It can be 
seen that 6 and d increase monotonically with |»?|. Neglecting terms involving 
negative powers of y?, the coefficients in the expression for S (v2) are given by 


2d d(v? + B, 6%) (v? — B, d2) = 283 _ 
; Psa) a ake 


d? (v2 + B, 62)2 — 62 (v2 — B, a2)? = (B:— Bs)? yo (B23 +28; Ba) (Bi—Bs) 28; B3_ v* 


ct yt Ba C2 y2 — Bg 


% 
C 


B, BR By B3 pt 
_, BaBS—{B3— (Bi —Bs) Ba} {Bs (B3—Br Ba) +83 Ba} »* 
63 (B, — Bs)? ca 


Thus, for large |»? 


, both coefficients in the expression for S(»2) are negative. 
Therefore it is certainly possible to choose a constant b,<0 such that for v?<},, 
the coefficients are negative and d, 5 increase without limit as y2 decreases without 
limit. Certainly there exists a constant b,< 6, such that for v?< by, cosh 61>2/)/2. 
There exists 6;<b, and an integer » such that for ve=bs, dl=(2n+%)x. At 
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once, S (b;)>0. Finally, there exists b,< 6; such that at »?=6,, dJ=[2(n+1)+2]a 
but S(b,)<0. Therefore there must be a solution v? to (9.6) and (9.12) in the 
range b,<y?<6,;<0. That is, the membrane is unstable for p,;<0. 

It remains to consider the region B,>0. From equation (9.6), we see that 
the nature of 6, and. 6, depends on 4,, A,, and v®. Suppose f, 6, —f3=0. If 
there is a solution »?<0 in this range, 6, and 6, must be of the following types: 

(7) 0: = 2, 

(77) 6, and 6, real, 

(271) 6, and 6, complex. 

For case (¢), equations (9.11) do not apply. Writing a general solution for 
case (7), we apply the boundary conditions and find that the resulting trans- 
cendental equation can have no solutions for which y? is negative. 

Solving (9.13) for v?, we obtain 


Diet lg at) (0, + 62) sinh 3 (d,— 5.) [= (0,— 6) sinh 3 (4, + 92) J 
ete Pes (6, + 62) sinh $(d,— 6.) 7+ (6,— 6,) sinh 3 (0,+ dy) 2° sda 


Suppose there is a solution y?<0 for which 6, and 6, are real, case (77). Both 
6, and 6, may be considered positive without loss in generality. Setting 


(0, — 02) 
se 


po PE ee, 


B, and By, are real, B,> 6,, and without loss in generality we may assume 
that B,>0. Then equation (9.14) becomes 


B, sinh ByF B,sinh By _ 


ee doit Fai) 
i Pr &1 Og B, sinh B, +B, sinh B, Pr Onis” 
where we have written 
p=) a B, B, _ [sinh B, sinh B, 
ns TH B,sinh B,+B,sinhB, | B, By 


since B, > B,>0, 
sinh B, sinh B 


2 <0 
B, By > 


and 7;*)<0. Therefore, since f, 6, 62> 0, v? is positive. This is a contradiction 
so that for B,>0 and f, 8, —f3> 0 there are no negative solutions for y? where 
6, and 6, are both real. 

Suppose »2< 0 and 6,, 6, are complex, case (777). Since 6,, 6, are then complex 
conjugates, we may write 


Gee Gnd Dawe Oo ealsedba> 


where }, and 6, may be taken to be positive without loss in generality. Sub- 
stituting into (9.14), we obtain 


ps 2, 72) 0, sin b, 1+ b, sinh by if 42 pe) TP). 
We By (0 b3) b, sin 6, / +b, sinh by 1 Bi ( 1 2) 141 
Now, 


= + {b;1| sin byl] — 6,7 sinh b, J} < Pal {01 — sinh 0,1} <0, 
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since 0<b,l<sinhd,/, b,1> | sin byl]. 


Therefore b, sin byl +b, sinh b,l> 0, 


b, sin b,/ — by sinh bl <0, 


which implies 7;;)<0. Hence v?>0, a contradiction. Thus the region f,> 0, 
B, Bs —B2> 0 is stable. 
Finally by a similar but more complicated analysis we can show that for 


B3> 0, 0< #3} BiBi< 4B, Pon*(F) , 


there are no negative solutions for y? but that for 


Bs>0, fP2—B,Bs> 4P,Bs72 (5) , 


there are solutions for »? which are negative. The details of the work will not 
be given here, as the method does not appear to be of general interest. 

In this way we are able to establish that the region defined by (7.10) is the 
stable region. 


10. Determination of the Natural Frequencies 


Equations (9.6) and (9.12) or (9.13) determine the natural frequencies of 
cylindrically symmetric vibrations, involving radial and axial displacements uw 
and w, for a circular cylindrical membrane of incompressible material which has 
been finitely extended and inflated. It does not seem possible to obtain an 
explicit expression for v? in terms of 4,, Ay, I/c, c even for the simplest form of 
the strain-energy function, that for a neo-Hookean material. To compute »? 
we may proceed in the following way. For particular extension ratios A,, A» 
and a particular length to radius ratio //c, we assume a value for vy? of the form 
t/c2 with t a positive constant. Substitution into (9.6) gives the values of 6, 
and 0, for this value of »*. Substitution of these values into (9.12) or (9.13) 
will tell whether the choice t/c? is close to a solution. Proceeding by trial and 
error we can obtain a sequence of values for y? which can be made to approach 
the true value to any degree of accuracy desired. From equation (6.3) we then 
obtain the natural frequencies w for fixed ends. This is a cumbersome procedure 
and it would be quite useful to develop an approximate method for the deter- 
mination of the natural frequencies. 

For the higher frequency modes it is to be expected that the influence of the 
boundary conditions at the ends of the cylinder on the values of the frequencies 
is small. By relaxing the condition of fixed ends, the computations are greatly 


simplified. One approach is to assume that the modes of vibration are approxi- 


mated by 2 
u = B,sin™ ; O° cos ot, 


(10.1) 

Ww =1B, Cos =~ se cosa t, 
where m is an integer and where B,, B, are constants. In these modes the ends 
remain circular (w=0 at the ends) but the axial displacement w at the ends is 
non-zero, and in the literature (see [9] for example) they are referred to as 
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“breathing modes”. The boundary conditions are referred to as “freely sup- 
ported”. Assuming solutions of the form 


m me O 
u = Bycos —7— cos w t, 


m mo? 
L 
leads to the same values for the frequencies. 


Substituting (10.1) into the equations of motion, we find that in order for 
B,, B, to be non-zero, vy? must be given by 


2 2 
QC ig ae Ba + (Bi + Bs) eee ae 


+ ( at (B: + Bs) War) 4 TF (1 Ba — B2 + By Bs TEE) 


Bynere M==1, 2)... ¢ 


w= B,sin cos wt 


(10.2) 


The stability region for the cylinder with freely-supported ends is found to be 


P20, PR P,B.<B.Bs qian (10.3) 


This region is contained within the stable region for fixed ends. 


Consider the solutions of (9.12) for large v?. In particular, we are interested 
in the range »2c?>f,. For this range, the quantities 6,, 6, defined in (9.10) 
are purely imaginary and (9.12) can be written S(v?)=0, where now 


S (v2) = [ei (v® — B, €3)? + €3(v? — B, e7)?] sin e, J sin e, 1 — 


(10.4) 
= 28 & (v2 — B, €2) (v2 — B, €3) [4 — cose, 1 cos eg/], 


and where 6,=76&, d2.=16&. We choose values for /,, A,,//c and assume a 
particular value for v?. Substitution into (9.6) gives e, and ¢,. Inserting »? and 
€, (v2), €(v?) into (10.4), we obtain S for this value of »2. Whenever S=0, this 
value for y? is a solution. 


Expanding ¢, and ¢, for large v?, we find that: 
(1) ¢, and e, both increase monotonically with v?, 
(2) 2, €. (v2 —B, €%) (v? — B, €5) is negative and of order 4, 


(3) ef (v? — B, e3)2-+ e3 (v2 — B, €7)? is positive and of order 7°. 


The first term in the expression for S vanishes whenever either e,/ or é5/ 
is equal to mz, m=1, 2,.... From (2) and (3) above, the sign of S, for large v?, 
can change only when ¢,/ or €,/ are close to ma, since except when sin ¢,/ sin €,/-~ 0, 
the first term in the expression for S predominates. From (9.6), when e¢,/ or 
&,l 1s equal to ma we have 


TF i. sae [Bi Ba Be (B, + Bs) cy? + 
{(B3 — Br Ba + (Bi + Bs) © »?)? — 4B, Bs (2% — By) 29738]. 


Solving this equation for y?, we again obtain equation (10.2). It follows, there- 
fore, that the frequencies for the breathing modes are good approximations to 
the frequencies for fixed ends when the frequencies are large. 
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Comparing (9.5) with (10.2), we find that on B,;=0, the frequencies are the 
same for both freely-supported and fixed ends. The same is true on the line 
B,=0. We thus expect that the freely-supported ends approximation will be 
useful in the region between f,=0 and 6,0 for small y? also. 


11. Numerical Results 


Fig. 6 compares the axially symmetric vibration frequencies as a function 
of A, for fixed and freely-supported end conditions when [’=0, J/c=6, and 
4,=1 (neo-Hookean material, no axial extension). We see that the frequencies 
coincide when £,;=0 and when f,=0 
as mentioned previously in Section 10. 
For the situation represented in Fig. 6, 
the frequencies for the two end con- 
ditions coincide at A, =1 and A, =1.315. 
Thus we expect that the agreement for 
intermediate values should be good even 
for small frequencies, as can be seen from 
Fig. 6 for the particular case A, = 1.25. 
Finally we notice that at 2,= 1.88, which 
from Fig.5 is where the tube becomes 
unstable, the lowest frequency goes to 
zero, as it should do. As is to be ex- 
pected, the agreement between the 
lower frequencies is poor for values of 
A, near the point where the cylinder 
becomes unstable. 


a, Po aa a eva as Figs. 7, 8, and 9 and the Table 
Fig. 6. Comparison of lowest frequencies for fixed and freely- indicate the changes in the lowest 
supported ends; //c=6, A4,=1, T=0. ——— Freely supported frequencies for freely- 
ends, — — fixed ends, x isolated values for fixed ends q £ vasre supported ends 


as the parameters are varied. 
Fig. 7 illustrates the variation of @c?@?/2C, with I. The two lowest frequencies 
are plotted against J’ for 4; =1.25, A,=1, and J/-=6. It can be seen that for 
fixed A,, A,, and l/c, the frequency parameter w? increases almost linearly with [. 


Table. Lowey Frequencies, Freely-Supported Ends; 4,=1.25, 4,=1, llc=2, P=0, 0.1 


9 ct w;)*/2C, 9 hat 2/2C, 
el r=o | r=o1 | T=0 T=0A 
| | | 
| 

1 2.48 | 2.74 7.23 8.40 
2 5.16 | 5.84 | 28.79 33.31 
3 9.60 | 10.98 | = 6473 74.87 
4 13.82 | *4$8.49 115.0 | | 433.0 
5 23.82 27.40 | 179.8 207.8 
10 | $046 | fod. * > eas or SIR Ve34 6 


In Fig. 8 @c?w}/2C,, where w, is the lowest frequency, is plotted against A, 
for A,=1, l/c=2, and ['=0, 0.4...The instability point is at a higher value of 
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A, for ’=0.1 than for the neo-Hookean material, /’=0. From equation (10.2), 
there are, for each integer m, two values for y? and therefore for w2. These values 
are denoted by w\*)? and w{>)? and for I/c=2, Ay=1, 4;=1.25, the Table gives 
the values of 0 c? wl +)2/2 Gy forwm=4 ; 2, 2.94. 


be 4, 5, and 10 and for /=0-and. 0.1, a 
l/c=2 
IT =0,07 
20 Ze) ; 
I= 07 
76 
76 20 
74 
I70 
> 72 15 | 
S 
aX 
as 70 
S 
aS 
6 70 
6 
¥ OS Ae 
2 
Zi 0 
0 OS 70 ee, IE PD) 7 75 2.0 a, aon 
Fig. 7. Two lowest frequencies as a function of I’ for freely- Fig. 8. Lowest frequency for freely-supported ends 


supported ends; //c=2, A,=1.25, A,=1 


Fig. 9 shows the variation of the lowest frequency with /,, for A,=1 and 
l/c = 2, 6, 10 in the case of a neo-Hookean material, /’=0. 


7 72 74 6. °~ «8 Doge Mae 


Fig. 9. Lowest frequency for freely-supported ends and various l/c 


For [’"=0, l/c=6, Apy=1, and A,=1.5 the first and second eigenfunctions 
UD, W® and U®, W® were computed for fixed ends. It was found that WY 
is an even function of £=#%—J/2 with no zeros except at the ends while U™ 
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is odd with one interior zero, at the mid-point =0. For the second eigen- 
vector, U® is even with no zeros except at the ends while W®) ,which is odd, 
has a zero at the mid-point. In Fig. 10, these eigenfunctions have been plotted 
against €. They have been suitably scaled so that the shapes of U®, WY, U®, 
and W’) may be compared. 


| 
“la -3U/6 = 1/4 1/8 
6? 
Fig. 10. Eigenfunctions, fixed ends, for two lowest frequencies with =O, I/c=6, 2,=1, 4,=1.5. ~ Lowest frequency, 


—— second lowest frequency 
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Greensche Tensoren und asymptotische Gesetze 
der eletromagnetischen Hohlraumschwingungen 


CLaus MULLER und Horst NIEMEYER 
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Einleitung 


Es sei G ein Gebiet des dreidimensionalen euklidischen Raumes, das mit 
homogener und isotroper Materie erfiillt ist und von einer ideal leitenden Rand- 
flache F begrenzt wird. Die Materialeigenschaften seien durch die Dielektrizitats- 
konstante « und die Permeabilitat w festgelegt. Sind ¢ und yw positiv reell, so 
existieren abzahlbar unendlich viele elektromagnetische Eigenschwingungen ©,,, 
H, und Eigenfrequenzen w,, die in G den (zeitunabhangigen) Maxwellschen 
Gleichungen 

Vx, +t, €¢€,=0 


Vx, iy TO, 9, = a) 


und auf F der Randbedingung nx ©,—p geniigen. Dabei ist n der ins AuBere 
von G gerichtete Normalenvektor auf FP. 

Die in dieser Arbeit behandelten Fragen traten erstmalig im Zusammenhang 
mit der Planckschen Strahlungsformel auf, bei deren Ableitung angenommen wird, 
daB sich die Eigenschwingungen fiir hohe Frequenzen im Mittel in allen Punkten 
gleichartig verhalten. Ein Ziel der vorliegenden Arbeit ist es, diese Annahmen 
zu prazisieren und im Rahmen der mathematischen Theorie elektromagnetischer 
Schwingungen streng zu beweisen. 

Bei den Randbedingungen beschranken wir uns auf 1x €=op, da bei physi- 
kalischen Uberlegungen diese Forderung bisher vorwiegend betrachtet wurde. 
Man kann jedoch erwarten, daB andere Randbedingungen, die heute gelegentlich 
in der Theorie der Hohlraumresonatoren auftreten, nach einem analogen Ver- 
fahren erfaBt werden k6nnen. 
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Fir die Eigenfrequenzen w, ergibt sich das zuerst von H. Wevt [17], [18] 
aufgestellte asymptotische Gesetz* 


wobei V das Volumen des betrachteten Gebietes G ist. 
Wir werden in dieser Arbeit iiber das Weylsche Gesetz hinaus fiir die durch 


lu 
G G 


normierten Eigenschwingungen die asymptotischen Relationen 


x |€,.@))? ~x (en)! a 


0<on<t 


Dd |B, (2)? ~ ele w)? 


0<an<t 


3 


Ba" 


(x €G, too) ableiten, aus denen sich das Weylsche Gesetz formal durch Inte- 
gration tiber G ergibt. Diese Gesetze kénnen als Formulierung der asymptotischen 
Ortsunabhangigkeit der Eigenschwingungen angesehen werden. 

Fiir die skalare Helmholtzsche Schwingungsgleichung wurden analoge Gesetz- 
maBigkeiten von T. CARLEMAN [1], [2] mit Hilfe Tauberscher Satze von HARDY 
und Litrr—Ewoop bewiesen. Zur Anwendung der Carlemanschen Methode auf 
elektromagnetische Schwingungen ist es notwendig, Greensche Tensoren zu den 
Maxwellschen Gleichungen zu konstruieren und ihr Verhalten fiir groBe imaginare 
Frequenzen zu untersuchen. Hierbei treten charakteristische Schwierigkeiten auf, 
da die Maxwellschen Gleichungen keine Grundlésungen besitzen, die wie 1/7 
singular werden. Die als Grundlésungen anzusehenden Dipolfelder 


E(t) =axV,@ 
SQ) =<meDat <P, (VY, ®) 


‘ ; kr 
mit einem konstanten Vektor a und ®— % 


: ,M=wre uw, OSarg&<a,r=|t—| 
werden vielmehr wie 1/7? bzw. 1/7 singular. 

Die Ubertragung der Idee der Greenschen Funktion auf Probleme elektro- 
magnetischer Schwingungen ist vielfach diskutiert worden. Der hier gegebene 
Nachweis der Existenz von Greenschen Tensoren zu den Maxwellschen Glei- 
chungen bezieht sich auf die Lésung des Randwertproblems, die der eine von 
uns 1952 gegeben hat [11], [12]. 

Wir beginnen daher in §§2 und 3 mit einer Untersuchung der Greenschen 
Tensoren 


G9) (4, k= 1,2), 


1 > 1 bedeutet die Anzahl der Eigenfrequenzen zwischen 0 und ¢. 


0<an<t 
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d.h. dreireihigen Matrizen, deren Spaltenvektoren die homogenen Maxwellschen 
Gleichungen und die Randbedingung erfiillen. Mit ihrer Hilfe laBt sich auch das 
inhomogene System 


VxG+¢oeE=j 
VxE—¢ouH=—j' nG 


1x G—=p auf F 
lé6sen, und zwar in der Form 


6= it f Oui+Greirar 


34WE 
G 
tistiat een ee 5 
D stall tag f Grit Geri aV 
G 


Hierbei sind die Integrale als Cauchysche Hauptwerte und die Produkte ,, j 
usw. im Sinne des Matrizenkalkiils zu verstehen. 

Um den Anschlu8 an die Carlemansche Theorie zu gewinnen, behandeln wir 
auch den Fall #=0. Obwohl das System der Maxwellschen Gleichungen fiir 
m=0 zerfallt, gelingt es, durch stetige Erginzung Greensche Tensoren in einem 
erweiterten Sinn zu konstruieren. Setzen wir e=u=1 und bezeichnen mit 
®&,,(x, 9; 27) die fir #=¢t (t=0) konstruierten Greenschen Tensoren, so ergeben 
sich fiir die Eigenschwingungen der Maxwellschen Gleichungen die Integral- 
gleichungen 

(wo, + 7?) iF Gait, ys ¢t) E(x) dV, = 40 4, 9H, (9) 


(w* aaako) 11 Bal (na); ¢T) 9, (t) dV, =42 40, g, (8) Ne 


Mit ihrer Hilfe laBt sich die Existenz von abzahlbar unendlich vielen Eigen- 
schwingungen beweisen. 


Der Grundgedanke des Carlemanschen Verfahrens ist jedoch erst anwendbar, 
wenn wir von dem iterierten Tensor 


1 wo 
Gr (t, 95 27) = gor | Silos 0) Ge, (3, 9; 2) AV, 
G 
ausgehen. Die Tensoren @,, und ,, werden wie 1/7? singular, der iterierte 


Tensor %, jedoch nur wie 1/r. Fiir ihn gilt die se ai tat vanes 


Ey (x 


wz am = 


Gz (x, 9; 0) — Ge (x Hip BAB j=2 > 


won >0 
aus der die asymptotischen Relationen nach dem Verfahren von T. CARLEMAN 
folgen. 

Bevor diese Uberlegungen in §§5 und 6 ausgefiihrt werden, ist es jedoch 
erforderlich, Abschatzungen fiir die kompensierenden Tensoren zu gewinnen. 

Die in §4 vorgenommenen Abschatzungen zeigen, daB die kompensierenden 
Tensoren nur dann singular werden, wenn Aufpunkt und Quellpunkt gegen den- 
selben Randpunkt von F streben. Analytisch laBt sich dies nach H. WeEyr [17], 


2A 
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[19] durch den Begriff des ,Lichtweges“ zwischen Aufpunkt und Quellpunkt 
fassen, also durch die Entfernung, die ein Lichtstrahl bei einmaliger Reflektion 
am Rand zuriicklegt. Die kompensierenden Tensoren lassen sich durch Potenzen 
des Lichtweges zwischen Aufpunkt und Quellpunkt abschatzen. 


§ 1. Problemstellung, Bezeichnungen 


Wir verwenden im folgenden die Bezeichnungsweisen der Vektoranalysis und 
der Matrizenrechnung nebeneinander. Dabei wollen wir uns Vektoren als Spalten- 
vektoren geschrieben denken. Das Transponieren von Matrizen und Vektoren 
deuten wir durch ~ an. 

Das Skalarprodukt a~ zweier Vektoren a und 6 werden wir jedoch auch 
mit a bezeichnen, wenn keine MiBverstandnisse auftreten. Unter dem Spalten- 
vektor a Xb verstehen wir das Vektorprodukt der Vektoren a und b, entsprechend 
unter Vxa die Rotation von a, usw. 

Fiir die in dieser Arbeit benutzten Begriffe Divergenz und Rotation von 
Vektorfeldern legen wir nicht die iiblichen differentiellen, sondern erweiterte, 
durch die physikalische Anschauung nahegelegte Definitionen zugrunde [/7], [13]. 


Definition 1. Das Feld v(r) set in der Umgebung des Punktes Yo stetig. Es 
sei G, eine gegen Yo konvergente Folge regulirer Gebiete! mit dem Volumen |G,]|| 
und der Randflache F,. Es set ferner nu der ins Aufere gerichtete Normalenvektor 
auf F. Existiert dann der Grenzwert 


11 lim Tay abyaF 
( ) Gy Xo IG, Ih ( ) 

Fy 
bzw. 
Ae lim rent nxv) dF 
( ) Gy Xo Gull. ( ) 


unabhingig von der Wahl der Folge G,, so stellt im Punkte Yq der Grenzwert (1) Vv 
der Grenzwert (2) Vxv dar. 


’ 


Fir stetig differenzierbare Felder stimmen diese Definitionen mit den tiblichen 
differentiellen iiberein. Man kann jedoch zeigen, daB es stetige, aber nicht 
differenzierbare Felder gibt, fiir die eine Divergenz im iiblichen Sinne nicht 
existiert, sondern nur eine Divergenz im Sinne der Definition 1. Da die bekannten 
Regeln der Vektoranalysis auch fiir die verallgemeinerten Begriffe Rotation und 
Divergenz gelten, ist es vorteilhaft, die Klasse der zugelassenen Felder zu er- 
weitern und solche stetigen Felder in die Betrachtungen einzubeziehen, fiir die 
Divergenz oder Rotation nur im erweiterten Sinn vorhanden sind. 


Es sei nun G ein endliches, regulares Gebiet, das von der orientierbaren, 
zusammenhangenden, dreimal stetig differenzierbaren Randflache F berandet 
wird. In G seien die Vektorfelder j, j’ und ihre Divergenzen Vj, Vj’ definiert 


und stetig. Entsprechend dem in der Einleitung skizzierten Weg behandeln wir 
zunachst folgendes Problem: 


Zum Begriff regulares Gebiet: [11], S. 20ff. Eine Folge regularer Gebiete G, mit 
%€G, und gegen Null strebenden Durchmessern ,konvergiert gegen ro‘. 
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Gesucht sind Felder ©, mit stetiger Rotation, die in G die Maxwellschen 
Gleichungen 


(1.3) Vx$+éweE=j 


VxE—¢ouH=—j’ 


mit komplexen Konstanten e, wu erfiillen und auf F die Randbedingung n x © =p 
befriedigen. 


Die Konstanten ¢ und yw schreiben wir in der Form 


éo" 


Oo 
es + =, {= by + i 
Dabei ist ¢) die Dielektrizitatskonstante, uw) die Permeabilitaét, o die elektrische 
Leitfahigkeit, o’ die magnetische Leitfahigkeit des mit homogener und isotroper 
Materie erfiillten Gebietes G. Die Konstanten é9, #9, ¢ und o’ sind reell und es gilt 


fg ea Wilp 10) gO 2 Oe OF 


In den physikalisch realisierten Fallen ist o’ =0. Die Frequenz w schranken wir 
durch die Forderung Inw=0 ein. Die Gln. (1.3) gehen aus den zeitabhangigen 
Maxwellschen Gleichungen hervor, wenn alle Felder den zeitabhangigen Faktor 
e*°* enthalten. 


Wir wollen diese Aufgabe mit Hilfe Greenscher Tensoren ;, (, )) in der Form 


CQ) = Ot f Gre iO)+ Gee v) i] 4% 


34WE 


(1.4) 


90) =O + ZF [Gul iW) + Gee iP O)]4% 
G 


lésen. Wie wir schon einleitend betont haben, kénnen wir wegen der Singulari- 
taten der Greenschen Tensoren die Integration nicht im gewohnlichen Sinn durch- 
fiihren, sondern nur im Sinne des Cauchyschen Hauptwertes eines bestimmten 
Integrals. Wir bendtigen daher 

Definition 2. Ist die Funktion f(x, y) wm Gebtet G mit Ausnahme des Punktes 
r= (y CG) tiberall stetig, so ist 


(1.5) Fit.y)dV¥,=lim f fe») dv, 
G 1" G-Ky 
falls dieser Grenzwert existiert. Dabet bezeichnet Ky eine ganz in G gelegene Kugel 
vom Radius r um den Punkt v. 
Ist f(x,y) (uneigentlich) integrierbar in G, z.B. |f(x, »)|=O(|z—»|~*) mit 


“<3, so ist fre b) v= J 10 y) aV,. 


§ 2. Singularitaten und Symmetrien Greenscher Tensoren 
Wir gehen von einer Integraldarstellung der durch Volumenstréme erzeugten 
elektromagnetischen Schwingungen aus, die uns die Greenschen Tensoren fiir 
den freien Raum liefert, und benutzen [/7] (Lemma 49 dort). 
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Lemma 1. Im reguldren Gebiet G seien j,Vj,j’, Vj’ stetig. Mit Vi=¢oe, 
V j'=¢wo' und 


tk 


P(e, y) = 2, r=|r—y|, A=aweu, OSarge<n 
Y 


werde das Feld 


Ser) | aN Na sen al ee 
(2.1) 


ee 


(i'n) Vy dF, 


(4 


5* (1) heh D4jxV) P+ £0, av 
G 


gebildet. 
Dann gilt in den inneren Punkten von G 


Vx$*¥+2meC®=j, VxeF—Zopnht=— yj, 
wahrend in den duperen Punkten 
VxX9* + soe =p, Vx*—Zouyp*—o 
erfillt wird. Ist auBerdem auf der Randfldache F 
(nj) = (nj) =0, 


so sind &* und $* in allen Punkten stetig. 


Wir schreiben die Felder Gl. (2.1) in eine unseren Zwecken besser angepaBte 
Form und beweisen: 


Lemma 2. Es ist 


Cot) ee eT fa f OHIO + ORO] A, 


o*t) = 371 react uatgencere 
G 


320 bu 


wober die Tensoren Wi, (x, y) durch die Forderung bestimmt sind, daB fiir jeden 
konstanten Vektor a 


GH (et, b)a=couPa+ — <. (a V, ®) 
Oz (rt, ))a=V,@8 xa 
(2.3) 31 (z, 9) 
Gf (rz, y)a=—V,®xa 
O45 (t, 9) a=¢ éoePa+% V,(aV,) 
ist. 


Man erkennt hier, daB die beiden ersten Ausdriicke das elektrische bzw. 
magnetische Feld eines infinitesimalen elektrischen Dipols vom Moment a an der 
Stelle ) darstellen, die beiden letzten Ausdriicke entsprechend das elektrische 
bzw. magnetische Feld eines magnetischen Dipols. 
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Zum Beweis von Lemma 2 betrachten wir 


[ @%%a%,= =P, [ i) oar, 
G 


(2.4) : 
=5%.| Poioar+ firs) an), 
F G 
Mithin ist 
Qs) f2noa%,+* fmnoan=7, f i%aa",, 
F G 


und es bleibt lediglich die Relation 


(2.6) hf (i 10) av ie) + fh i% ®) ay 
G 


nachzuweisen'. Dazu legen wir um den inneren Punkt xy)CG eine Kugel Kj, 


vom Radius t, die noch ganz in G enthalten ist, und bilden mit |r —r)|< = 


ef iV; ®) AV, =V, S iV, DAV, +V, SUR) a 
G 


é- Ky, Ky, 


0 


= JS Veli bie) VO) dN hike dV, 


(2.7) = /VUGV,®)dV,—V, f (ni) Oar, +¥, f Vi) av, 
G—Ky, In tol =7 Ky, 
= JVM Pav + f nih Oak, +f Vi) Vay, 
(GS —Ky, [y— Uol = =t Kr, 


Setzen wir nunmehr y=yz, und lassen t->0 streben, so folgt wegen 


; AH +), 
him (uj) KP a, = — — j 0) 
Ip—tol=2 
und 


S Vj) V,@ aV,| = 0 (z) 


Ke 


die Behauptung Gl. (2.6). 

Wir kénnen nach Lemma 2 die Tensoren ©}, als Greensche Tensoren fiir den 
ganzen Raum auffassen. Mit ihrer Hilfe konnen wir nun auch die Greenschen 
Tensoren eines endlichen Gebietes definieren. 

Definition 3. Es sei w keine Eigenfrequenz des Gebietes G. Die Tensoren 

in (t, 9) heiBen Greensche Tensoren des Gebietes G, wenn folgende Eigenschajten 
erfiillt sind : 

1. [G;,(x, 9) — GA (e, y)] st in jeder Komponente als Funktion von x stetig 
differenzierbar in G. 


110]; S235 it. 
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_ Mit jedem konstanten Vektor a erfiillen die Felder &,,a und &5,a als 
Pacer von x inG fiir y+-) die homogenen Maxwellschen Gleichungen, d.h. 


V, x Gy ,0 a PORT AOP TE =D 
V, x G2,a+ 40 E G,,a= D. 
3. Es gelten die Randbedingungen 


xX (Gz, 0) =0 
2) ie are es OP) 

Bevor wir in §3 die Existenz und Eindeutigkeit der Greenschen Tensoren 
beweisen, wollen wir einige allgemeine Eigenschaften der Greenschen Tensoren 
herausstellen. Wir werden diese Eigenschatten bei dem Existenznachweis nicht 
bendotigen. 


Satz 1. Die Greenschen Tensoren erfiillen die Symmetriebeziehungen 


Oni (r, 9) = G11 (0, z) 
Gro (v, 9) =— Gor (d, 2) 
G22 (tz, 9) = Gee (b, z)- 


Zum Beweis von Satz 1 wollen wir zunachst zwei Hilfssdtze herleiten. Es 
sei Fy die Flache einer Kugel von Radiusy um den Punkt y, n der ins Innere 
gerichteten Normalenvektor, a ein beliebiger konstanter Vektor und f(r) ein 
stetiges Vektorfeld bzw. ein Feld mit stetiger Rotation. Dann gilt: 


Lemma 3. 
lim J w[f(t) x(%® xa)] dF, =~ af(y). 
B 
Lemma 4. 
lim { n[f() )xV,x(V,® xa)] dF = =* a(Vx#) (0). 
Fy 


Beweis zu Lemma 3. Es ist (’ = d/d7) 


(2.8) Lo ("\'n, 

folglich 

(29) f uff) xV,Oxa)]aR= ("2") f [ova) (nf) —aflan, 
By BS 


und wir erhalten 


(2.10) lim (yf (af) dF. = —4zxaF(y), 


(2.11) Him (2°) f (v0) (nf) a = — lim 4 (na) (nf) dF, = — 4” af(y). 
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Beweis zu Lemma 4. 


lim [ 1G) x, x(%,®xa)] dF, = lim [ a0F@) xP,(a V,.®)| dF, 


Fy Fy 
’ F 4 ZAP \e 
(2.12) = lim © (2) | (f(z) adF, 
Fy 


Zum. Beweis von Satz 1 bilden wir nunmehr mit konstanten Vektoren a, b 
Gi2.—= Gi (t,9)0, YEG, HEF, 

Gin = Win leah RPSEG, FEF, 425b, 

und betrachtet den Ausdruck 


(2.13) 


(2.14) i [9i4°V_ XV Xin — Sim Ve XV XG;n] 2V,- 
, G—Ky—K3 
Da der Integrand wegen 
(2.15) V, XV, XG¢4 = OEM Gin =F Gir 
Ve XV eX Ope =O" E 1 Gin — FA? Wie 
verschwindet, erhalten wir nach dem GauBschen Satz 
(2.16) =e J 1 [Gi'm XV Gin — Gir XV_ X Gin] 4. 
F+Rj + 


Setzen wir 7=/=1, so folgt aus der Randbedingung 


(2.17) it Sp) eet pee oa 
(2.18) UF J natn Ve X G12 — Gan XV_X Gin] af,. 
By+ Fy 


Fir k=m=1 ergibt sich als Grenzwert fiir 0 nach Lemma 3 und 4 


tim f n[oh x, xu — 911 x%X 0h] aR, 
Fy 


=4o p lim [ [ot x (VP xa) dF + 


a 
(2.19) jn ele nl(V.xgh) XV, x (VB xa)] dE, 
Fy 


pe Swe * é oA, ; 
page OP Ga) 5 a 


V, XV, x Ort 
= 4270 (agi), 
ferner entsprechend 


(2.20) lim J ufgt XV, Xo. — 011*V, xg dR = — 4am pbgy. 
—- rR 
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Aus Gl. (18) folgt 


(2.21) a~ 11 (0, 3) 6 = b~ G1 (5, Y) a 
Gl. (21) gilt fiir beliebige Vektoren a, b, es ist daher 
(2.22) Ori (, 3) = G11 (5, 9)- 


Setzen wir in Gl. (18) k=m=2, so erhalten wir analog durch den Grenziibergang 
7—0 


(2.23) G22 (y, 3) = G20(5, y), 
und fiir k=1, m=2 folgt aus Gl. (18) schlieBlich 


Auf adhnliche Weise ]aBt sich folgendes Resultat gewinnen: 
Satz 2. Ist # reell oder rein imaginar, d.h. R=wre w=, so ist 
G1 0(t, ») = Gre (e, 9) 
Gait, 9) = Bai, 9) 
oe G11 (t, 9) =— @ € G1 (z, ) 
ot Gyo(t,) =— ou Gao(t, y), 


die Tensoren S15, Sg, sind tn diesem Fall also reell, die Tensoren we®,, und 
w@USs> rein tmaginar. 


Beweis. Es ist wegen Gl. (15) 
if [Sin- VxV.XGiim — Dim V.XVieX Gi] dV, 


25 G—Ky—K3 = a 
Cm =f 4) 1 area ee = os 
G—Ky—K; 


und fiir 7—/—1 folgt unter Beachtung der Randbedingungen 


(2.26) OS J [Gin XV_X 1% — Sin XV_ X Gi] dF,. 
Byt+ ky 


Beriicksichtigen wir nun, da8 analog zu Lemma 3 und 4 


(2.27) lim [ n[i@) xQ%® xa)] 4k = ie af(y), 
re 

(2.28) aa 1 [() xXYUx(V, @ xa) 3 eV xi) (y) 
Ry 


gilt, so erhalten wir zunachst im Fall k=1, m=2 


(2.29) lim Julgts xV, X11 — Sia XV XGf2] dF = 4a ¢o wage, 
Fy 
(2.30) iim Jn [gts xV,X 911 — O11 XV, X 9%] dF, =42b(V, Xo), 


3 
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d.h. es gilt 
G12 (9, 3) = — Gai 5, 0) = Gy 2 (y, 3) 
Goi (9, 5) = Ger (v9, 3). 


Entsprechend erhalten wir fiir k=m=1 und k=m=2 die beiden anderen Be- 
hauptungen. 


(2.31) 


§ 3. Existenz der Greenschen Tensoren 
Wir wollen in diesem Abschnitt die Existenz der Greenschen Tensoren mit 
den in §2 geforderten Eigenschaften nachweisen, falls m keine Eigenfrequenz 
unseres Randwertproblems ist, d.h. falls es keine von Null verschiedenen und im 
abgeschlossenen Bereich G stetige Lésungen (,, Eigenschwingungen“) ©, des 
homogenen Systems Aaa SRS 


Vx$+¢meC=o 
t= pat 
gibt. 
Wir zeigen zunachst ([/1], S. 319) 
Satz 3. Eigenfrequenzen treten nur dann auf, wenn « und u reell sind. Die 
Eigenfrequenzen sind ebenfalls reell. 
Zum. Beweis nehmen wir an, ©, § seien Eigenschwingungen und bilden 


(3.1) JUEXS) i fm Ex) dF. 

Es folgt 

(3.2) =f (@eEE—apnh $) dV, 
G 


und mit m=w,+¢, durch Zerlegung in Real- und Imaginarteil 
G3) J(u 8 9)dV =f (oEGAV, 
(3.4) O= f [(o’ +0; Mo) 8 H + (6 + oy &) EE] AV. 

G 


Die GréBen o, 0’, @1, &, Mo hatten wir als nicht negativ angenommen. Gl. (3.4) 
kann daher nur bestehen, talls c=o’=@,=0 ist. Damit ist Satz 3 bewiesen. 
Wir setzen nunmehr voraus, w sei keine Eigenfrequenz, und wollen die 


Exist der Tensoren 
xistenz ®;,(t, 9) = Gj, (e, 9) — GA 9) 


nachweisen, die als die kompensierenden Anteile der Greenschen Tensoren nach 
Definition 3 folgende Eigenschaften besitzen: 

1. ®;,(x, )) ist in jeder Komponente als Funktion von ¢ stetig differenzierbar 
in G. 

2. Die Felder ®,;(r, y)a und &,, a erfiillen als Funktion von y tiberall in G 
die homogenen Maxwellschen Gleichungen, d.h. 


VX 81,4 — 4 Ww Xy,a=0 


6-5) VX Rg,0 + ¢o € Ry, A= (Reyf2). 
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3. Es gelten auf F die Randbedingungen 
(3.6) nx ®, a =—n(t') x | do Ba+ ——V, (aV, )), 


(3.7) nXR,a=n(t’) X Vy B xa], 
CEP; yoGy Soer. 


Die Frage nach der Existenz der Tensoren 8;,(r,)) ist mithin ein typisches 
Innenraumproblem, das wir allgemein so formulieren kénnen: 

Gesucht sind Felder ©, §, die in G stetig sind, stetige Rotation besitzen, den 
ee eee Vx$+éweE=0, VxE—¢ouH=o 
gentigen und auf F die Randbedingung 
(3.8) nrxE=q 


erfiillen, wobei g ein vorgegebenes, differenzierbares Flachenfeld auf F ist, dessen 
Ableitungen hdélderstetig sind. Die Lésung des Innenraumproblems und des ent- 
sprechenden AuSenraumproblems ist in einer Arbeit des einen von uns ({/2], [11], 
S. 317ff.), gegeben worden. Wir wollen hier die Ergebnisse insoweit beschreiben, 
als wir sie spater noch benétigen werden. Die Lésungsfelder ©, 5 kénnen danach 
durch ein Flachenfeld t auf F erzeugt werden, das ebenfalls hdlderstetige Ab- 
leitungen und damit auch eine hélderstetige Flachendivergenz V, {= @, besitzt. 
Die Felder 


Eq) =— | (F0)x%%) ak 
(3.9) arr ; 
ce [so eDj + 2.7, 0) ak, 


erfiillen unter den genannten Voraussetzungen iiber f fiir alle r¢ F die homogenen 
Maxwellschen Gleichungen, auBerdem existieren bei Annaherung an die Flache F 
die Grenzwerte © und § , und es gilt (x’ € F) 


(3.10) nx €W)=Fie) — ge fre) x (FO) x%) ah, 


cat) 
F 


Fiir das erzeugende Feld f erhalten wir also die Integralgleichung 


(3.11) Fe) — 5 | mle’) x (Fy) x Hy ®) dF, = 2900’) 
symbolisch ‘ 
(E — M){=2g. 
Fithren wir als Norm des Elementes f 
(3.12) [7] =supr |f@’)| 
ver 


ein, so ist die Transformation M im Banach-Raum aller stetigen Flachenfelder 
auf £* vollstetig. Gl. (14) ist daher fiir alle g eindeutig lésbar, falls die homogene 
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Gleichung 
(3.13) (E — M){,=0 


keine nicht-trivialen Lésungen besitzt. Dies ist genau dann der Fall, wenn 
keine Eigenfrequenz ist. Unter den angegebenen Voraussetzungen iiber F und 
g sind die Lésungen der Gl. (11) hdlderstetig differenzierbar. Damit ist die An- 
nahme der Randwerte garantiert. Die Felder §t;,a sind auch eindeutig bestimmt, 
falls wm kein Eigenwert ist. Es ergibt sich also: 


Satz 4. Ist w keine Eigenfrequenz des Randwertproblems, so existieren die 
Tensoren &;;,(£,)) bzw. ©; , (x,y) mit den in Definition 3 angegebenen Eigenschaften. 
Sie sind eindeutig bestimmt. 

Wir wollen zum AbschluB dieses Paragraphen den bisher aus unseren Be- 
trachtungen ausgeschlossenen Fall w=0 diskutieren und zeigen, daB auch in 
diesem Fall Greensche Tensoren in einem erweiterten Sinn existieren. Formales 
Einsetzen von w=0 in Definition 3 ergibt keine sinnvollen Relationen, da schon 
die singuléren Anteile @7, (xz, y) und G¥g(r, y) der Tensoren ©,, und ,, fiir 
@=O nicht existieren. Wohl aber besitzen die Ausdriicke w®*, und wG#, fiir 
@->0 einen Sinn. Die definierenden Eigenschaften der ausgearteten Greenschen 
Tensoren fiir #=0, die wir mit &,,(z, ); 0) bezeichnen wollen, kénnen wir uns 
so erhalten denken, daB fiir w—>0 die Tensoren &,, und ,, gegen die ent- 
sprechenden ausgearteten Tensoren streben, wahrend bei den Tensoren ®&,, und 
®,. die Ausdriicke w@®,, und w®,,. gegen die entsprechenden ausgearteten 
Tensoren streben. Wir bilden daher 


Definition 4. Als ausgeartete Tensoren zur Frequenz w=0 bezeichnen wir 
Tensoren &;;,(x, 9; 0) mit folgenden Eigenschaften: 


1. Die Felder 
®yo(z, ); 0) a = Gyo (z, 9; 0) a+ V,— Xa 
®31(E,9; 0) a= Ga1(¢,950) a= Vy xa 
(3.14) 
z(t, 9; 0) a= G,, (, 9; 0) a — <P, [a V+ 
(k=1,2; r=|¢—»)) 
sind bet konstantem a stetig differenzierbar in G. 
2. Sie erfiillen in G die Differentialglerchungen 
VX Kya a = 7Ry oA, V, R14 = 0 
VX Xy2a = 0 
VX Ry, a= — ¢8,, 4, V,R2,a= 0 
VX 110 = 0. 


(3.15) 


3. Es gelten die Randbedingungen 


nx©,,a=D, n:-©,,0a=0 


16 
pre) nxG,,a=—09, auf Ff. 
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Wir wollen Existenz und Eindeutigkeit auch dieser ausgearteten Tensoren 


nachweisen und beweisen zunachst : 


Lemma 5. Die Integralgleichung i 
/ A . ii , 
GA) te = fae) x(to)xHl)am  (r=le’—9l) 
F 


besitzt fiir Im(A) +0 und fiir |A| <1 keine Eigenlisungen. 
Zum Beweis betrachten wir das mit einer Lésung f(r’) der Integralgleichung 


(17) gebildete Feld 
A 
(3.18) Ce) =A [ (to) x% 7) ah. 
F | 


Es gelten die ,,Sprungrelationen“ [//] 
U =f) , , /, , 
(3.19) nx€ (e')=+5“te), ax€ ¢') =H “te, 


wobei der Indexz Annadherung an F von innen, der Index a Annaherung an F 
von auBen bedeutet. Wegen 


Z 1 
Ce) = {Mx f ty) aK, 
F 
gilt ferner 
“ A : 
(3.20) A * 1 
F 


und wir erhalten folgende Relationen 


(3.21) 


Nunmehr bilden wir die Ausdriicke 


LV(EXVx® dv =f \VxE Pav 
i Gi 


=f(nx€)Vx€) dF, 
F >t 


=i 


(3.22) f\vxepav 454 [tS aP. 
x : >i 
Entsprechend erhalten wir 
(3.23) fivxePav = = a [ts dF. 
G; fares 


Hierbei haben wir zu beriicksichtigen, daB fiir 


le] = RS oo Cl Ol ral RRC = Olan) 
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ist. Zunachst erstrecken wir daher das Volumenintegral tiber den AuBenraum G, 
mit der Einschrankung |r|<R und erhalten (bei geniigend groBem R) zusitzlich 
ein Randintegral iiber die Kugelflache |r|=R. Dieses verschwindet aber fiir 
Ro wie O(1/R?). Nunmehr erhalten wir aus Gl. (22, 23) 


(3.24) (1+ 4) J|IVxEPPaVv + (1— a) f|7xEPav =o. 


Ist Im(A)=-0 oder |A|<1, so folgt unmittelbar Vx €=0 in G; und G,. Ist 


A=\ Tt , so verschwindet Vx © in ie 


Wir behandeln zunachst den ersten Fall (|A|<1 oder Im (A) ==0). Hier ergibt 
sich aus § =O sofort V)f=0 oder 


4 


a 


(3.25) fto) V1. dF,=0. 

F 
Setzen wir zur Abkiirzung 

ben A 1 
v= — | to) 48, 
F 

so ist 
(3.20) CSPxu— Vt 0, “A=; 


Bilden wir daher wiederum 


(3.27) fV(oxVxb) dV = f |Vxv|?adV = f (nxxd,) 0dF 
G; G; F 


f lVxvtav =— >" [tear, 
G; F 


und entsprechend 

(3.28) f\rxoPav=*5* [ toaF, 
Ga F 

so folgt analog 

(3.29) (1+ a) f|\Vxv|?dV + (1 — a) f|Vxv|24dV =0, 
G; Gg 

d.h. unter den obigen Voraussetzungen 

(3.30) Vxv=o inG, undG,. 


Nach Gl. (19) folgt f(x’) =o. 


Wir miissen noch die Falle A= +1 erledigen. Ist A= ee 


so verschwindet 


oe AG 
in lee Dann folgt aber wegen n x hee =p und V§$=VxH=0, daB  auchin ie 
Peieciy iden mu. Denn § laBt sich als Gradient einer Potentialfunktion an- 
setzen, die auf F konstante Randwerte annimmt. Wie oben kénnen wir jetzt 


schlieBen, daB Vx in Fe 


a 


? verschwindet. Mithin la8t sich » in ie als Gradient 


a 
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einer Potentialfunktion darstellen: »= VU. Das Feld v nimmt seine Randwerte 
stetig an, und wir erhalten 


(3.31) ffodF=f tr) Ve UdEy = [VU dF = — f Vt) UdF = 
F F F 


also ist auch VXb=0 in fee oder f(r’) =0, was zu beweisen war. 
Ti 


Lemma 5 setzt uns imstande, die Existenz ausgearteter Tensoren (fiir w =0) 


nachzuweisen. 


Satz 5. Es existieren Greensche Tensoren &;,(r, 9; 0) mit den in Definition 4 
geforderten Eigenschaften. Diese Tensoren sind eindeutig bestimmt. Es gelten die 


Symmetriebezxiehungen 
Sra (x, 9; 0) = — Ga (0, F; 0) 
(3.32) ia 
Grn (e950) =Gyay,t50) (k= 1, 2). 


Beweis. Wir beschaftigen uns zunachst mit den Feldern ©,,a und ,,a 
und erzeugen die kompensierenden Felder gemaB 


(3.33) Qea=-f(tG.o)xRt)\aR  (r=|z—s 
F 


% 


durch ein Flachenfeld f auf F, fiir das wir aus der Randbedingung n x G,.a=0 
die Integralgleichung 


(3.34) —2n(r’) x Dy X Ter = f(r’, y) — ay [ x(txQ >) ak, 
F 

erhalten. Die homogene Integralgleichung besitzt nach Lemma 5 keine nicht- 

trivialen Lésungen, die inhomogene Integralgleichung ist stets eindeutig lésbar, 

das Lésungsfeld f(x’, y), »€G, y¢F iitberdies hélderstetig differenzierbar. Die 

Tensoren ®,.(x, ); 0) und G,.(r, y; 0) sind auf diese Weise eindeutig bestimmt. 

Wir bemerken noch, da8 nach Gl. (33) 


jt 1 ta. 
(3.35) Ryo 0 VX me fe : dl, 
F 
ist. Es gibt also ein fiir y=) stetiges Vektorfeld tw(r, ), so daB 
(3.36) ®,.a=V,xiw(z, y) 
ist. 


Um die Existenz der ausgearteten Tensoren &,, und ©,, nachzuweisen, 
bestimmen wir zunachst ein Feld » mit folgenden Eigenschaften: 


(a) aires ist als Funktion von ¢ zweimal stetig differenzierbar in G, 
(b) fiiry +n ist Ap=o, 
(c) auf F gilt Vy>=0, nxb=p, 


wobei hier die Grenzwerte bei der Annaherung an F von innen her zu nehmen sind. 
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Das Feld v ist eindeutig bestimmt, denn die Differenz § zweier Lésungsfelder 
ware in G zweimal stetig differenzierbar und hatte die Eigenschaften (b) und (c). 
Mithin ware 

SV(9V5+5xVxo) dV = sd [(Vh)? + Vxh)?] dV 
(3-37) i i? 
sal [(1) Vo + (n xh) Vxh] dF =0, 


d.h. ¥h=0 und Vxh=p in G. Es wire also § der Gradient einer Potential- 
funktion, die wegen nxh=o auf F konstant sein miiBte. Daraus folgt aber 
§=0inG. 

Fir das mithin eindeutig bestimmte Feld ») machen wir mit einer noch unbe- 
kannten Belegung o (3, )) und einem noch unbekannten Flachenfeld j (3, y) (3 €F) 
den Ansatz: (r=|3—1]) 

(3.38) pee ieee! fen tah +s | (ixyt)an 
F 


| Bi 


der ersichtlich die Eigenschaften (a) und (b) besitzt. Wir bestimmen zunachst 
die Belegung 0 so, daB Vv auf F verschwindet. Es ist in G 


= 1 1 a 4 
(3.39) V8 =a Vie = bers dF, 
F 


Lassen wir den Punkt yx’ €F riicken, so ergibt sich nach den bekannten Sprung- 
relationen fiir das Dipolpotential fiir die Belegung o die Integralgleichung 


(3.40) aa f 695, FiR=— 0h, 


deren Lésbarkeit aus der Potentialtheorie wohlbekannt ist (Dirichletsches Pro- 
blem). Die Lésung 0 (r’, y) ist als Funktion von y’ (fiir y¢F) hélderstetig diffe- 
renzierbar [6], S. 212. Weiterhin bestimmen wir jetzt j so, daB auch die Rand- 
bedingung n xv=vp auf F erfiillt wird. Dies liefert uns fiir die Integralgleichung 


)— ax J 70) (i 6» 0) XV, 7 )a% 
=n x(a tag [end ren), 


deren Lésbarkeit durch Lemma 5 gesichert ist. Das Lésungsfeld j(z’, )) ist 
holderstetig differenzierbar. Damit ist die Existenz von » nachgewiesen. 


(3.44) 


Aus dem Feld » erhalten wir nun ©,, a und ,, a durch 
(3.42) Goya =V, x, 
(3.43) @,,0=i(V,xV, xb). 


Die Felder V, xv und i(V, xV, xv) besitzen naémlich zunachst die nach Gl. (3.14) 
geforderten Singularitaten. Ferner ist wegen 4,)=0: 


(3.44) V,xV; XV, x0 =V, x (VV, 0) — 4;o) = 09 
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also sind auch die Differentialgleichungen (15) befriedigt. Da die Felder V, xv 
und V,xV, xb ihre Randwerte gleichmaBig annehmen, gilt wegen n xv=0 aid 
V,, b=0 aut F schlieBlich n(V, xv) =0 und nx (V, XV, xv) =n xV,(V,0) =0 auf F. 
ite Gl. (37) liest man im iibrigen ab, daB auch die Felder &,, a und &,, a durch 
die Relationen Gl. (14, 15, 16) eindeutig bestimmt sind. 

Wir beweisen zum AbschluB die Symmetriebeziehungen und bilden mit kon- ~ 
stanten Vektoren a und b 


(3.45) Gin = Gi, (0, 0; O) a, Gip = Gi (x, 53 0) 6, 

sowie den Ausdruck 

(3.46) J V(gieXe1 — 4-0 Xgz2—4- w* Xg,,) aV,. 
G—Ky—K; 


Nach Gl. (33, 39, 45) ist 
V(gi2 921) =* 921 922 + 7 9f29i1, 
(3.47) — iV(v X93) = — 1 Go1 932, 
—1V(w* X O11) = — 7972 G11, 
der Integrand verschwindet also in G — Kj — Kj und nach dem GauBschen Satz 
folgt 


F+Rj+ KR 


Beriicksichtigt man, da neben den Randbedingungen Gl. (16) auch nxv=0 ist 
auf F, so erkennt man, daB das Integral itber F verschwindet. Die restlichen 
Integrale tiber A) und F” berechnen wir mit Hilfe von Lemma 3 und 4. Es ist 


lim | 1(gr2Xgo1 — 1 Xg22— 1 10* Xg,,) dF, 
i> : 


(3.49) = lim nats x (F; . xa) + w* x Vx (Fp xa)] dF, 
Fy 
= 42 a~ &,5(y, 3; 0) b 


und entsprechend 


lim Jn (gf XGo1— 70 X 22 — ¢ w* X qq) dF, 
Ober 
F 


” = 470 6~ Gi (5,9; 0) a 


Insgesamt ergibt sich daher 

(3.50) (1 9(0, 5; 0) = — G23, v; 0). 

Auf den analog vorzunehmenden Beweis der Symmetriebeziehung 
Gen (e, 9; 0) = Gy, (0,450) k=1,2, 


die wir weiterhin nicht benotigen, verzichten wir hier. 
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§ 4. Das Verhalten der kompensierenden Tensoren am Rand 
Wir wollen das Verhalten der kompensierenden Tensoren ®;, (x, y) am Rand 
untersuchen. Ist )€G, so sind sie fiir alle aus y dem abgeschlossenen Bereich G 
stetig. Wegen der Symmetriebeziehung (Satz 1) sind sie fiir r€G auch stetig 
fiir alle yCG. Als ein’Ergebnis dieses Paragraphen wird sich herausstellen, daB 
sie nur dann singular werden, wenn y und y gegen denselben Randpunkt auf F 


streben. Gema8 Gl. (3.9) werden die Felder 8;,(r, ))a durch Flachenfelder f, 
auf F erzeugt, die als Lésungen der Integralgleichung (3.11): 


(4.1) (E—M)tp=G, (k=1,2) 
gegeben sind. Fiir die Felder 

(4.2) ile > Ge oe MO oe fh=t 
erhalten wir die Integralgleichung 

(4.3) (E — M) i; = M" g, oe (UNG anol 


Dabei ist nach Gl. (3.6) und (3.7) 
Q(t’, 0) = — 2n(t’) x vom Da+ Ve (a Vy ®)| 
a Go(v’, ) = 2n(v’) x Vy ® xa). 


Wahrend die Flachenfelder g, fiir yz’ wie |x’ —|~° bzw. |x’ —y|-® singular 


werden, sind die Felder M@” q, fiir geniigend groBes n iiberall stetig, wie wir weiter 
unten zeigen werden. Folglich ist dann auch f% stetig, und die Singularitaten 


von f, kénnen durch 


(4.4) lee oe EG ee ee 
beschrieben werden, wobei gq, der am starksten singular werdende Anteil ist und 
die Singularitaten der Felder Mq,, M*q,, ... von immer geringerer Ordnung sind. 


Wir werden zunichst einige Hilfssatze itiber das Randverhalten transformierter 
Flachenfelder Mg beweisen. Dazu stellen wir das Flachenfeld Mg mit Hilfe einer 
dreireihigen Kernmatrix I(x’, ie in der Form 


(4.5) M g(t = J Mi (v’, 3) 93,0) dF, 


dar. In dem speziellen Fall der mae Gl. (3.6) definierten Integraltransformation 
ist 
(4.6) 22M(x’,3) =[n(z’) V, Gi, x’)] € — V, OG, v’) [{n(v’) xn )} xn) ]~. 


Hierbei ist & die dreireihige Einheitsmatrix. Im folgenden werden wir lediglich 
O-Abschatzungen geben. Daher wollen wir der Einfachheit halber vereinbaren, 
daB |%|=O(...) bedeuten soll, daB jedes Element a;, der Matrix & dieser Ab- 
schatzung gentgt. 

Wir wollen zunadchst Abschatzungen fiir das transformierte Flachenfeld 
M q(x’, )) gewinnen. Dazu zerlegen wir bei gegebenem r’€F undy€G die ae F 
in drei Teile F,, F und F,. Zum ersten (4) —3|\< 
(F,) alle Punkte 3¢€F mit |3—y|< zu 1, alle 
iibrigen Punkte von F. Die Punktmenge /, kann auch leer sein. 


i) 
bo 
* 
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Es gilt dann 

Lemma 6. Die Kernmatrix M(x’, 3) und das Vektorfeld 4 (3, y) seven definiert 
fir v' CF, 3€F, yEG und stetig fir x’ +3 bzw. 3+. Es seien folgende Voraus- 
setzungen erfuillt : 


(4) |oG.9|=O(lg—9l%), 0Se52, 


(2) | F960) aF,] 01) néG, 7'€F, 


(3) ene te sale 
(Aye te ay — es a fe a fiir |v al<sle’ b’|. 
Dann 1st 
OU) ane fe si 
(4.7) y Wee’, 3) 9 (3, 9) O(log =e falls ey 
O(1) Oe 


Beweis. Wir schreiben { = at + J + a und schatzen die Anteile einzeln ab. 
F 


=o. eariSr |a—n|* if) 


a 1 1 
(4.8) ae Ber iF) 


=O(|r’—5|~**") baw. O(1) 


Auf F ist |3—y|=4 


fee 3) 9(3, 9) a F; 


falls a >1 bzw. a<1 ist. 
Ist 3 €F3, so gilt 
Vey Se 


2 [2° he] 


und 
bby’ 
Gui 
( 7 = Pose = to = 


der Quotient lat" ist also auf F durch die positiven Zahlen = und 6 nach 


unten und oben beschrankt. Es ergibt sich daher 


saldlve ry lap] an) 
(4.10) =o{ f rpctpnth = of. [ e-sent) 


=0(|r’—y|-**1) baw. oft. a bzw. O(1) 


ine 3) 9(3, )) )aF, 


falls a>1 bzw. «n=1 bzw. a <1 ist. 
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Bei den Abschatzungen in Gl. (8) und (9) haben wir die Voraussetzungen (2) 
und (4) nicht verwendet. Lediglich bei der entsprechenden Behandlung der 
Integrale iiber /, bendtigen wir sie. 


Es sei y’ der FuBpunkt des Lotes von y auf F, y’ CK, der fiir geniigend kleine 
ae yr’ —y| eindeutig bestimmt ist. Fir 3€F, ist |r’ —3| 2 $|x’ —y|, speziell 
|x’ =+|r’—y|. Wir schreiben nun 


J Mx’, 3) 9%, v) 4, 
(4. 11) Bs Fr yY , / 
= ay) ) S96, ) ak, + J [M (x, 3) —Me’, olay) 4H, 
und erhalten 


(4.12) Pits 0) J 9G, 0) dF] =0(|x —y|-*#) baw. O(1), 


sowie, unter Benutzung von Voraussetzung (4), 


la—n’| 
sage [7 |a—n|* ir) 


aL aR) =O (5 2") bzw. OU). 


eae y)]9 (0) dF 


Wir diirfen SS ae 6 hier anwenden, da fiir zwei beliebige Punkte 3, €/, 
und 3,€F3: |3:—32|S3|a:—2’| gilt. Wir bemerken noch, daB Voraussetzung (2) 
automatisch erfiillt A falls «<2 ist. Setzen wir an Stelle von (2) 


| I a() iF —O(1) (fiir alle y€G) 
F 


voraus, so gelten die Aussagen unseres Lemmas auch fiir 0<a<}3, jedoch fiir 
a—=2 nur mit dem Ergebnis O(|z’ 1). Im iibrigen vgl. man 
H. Weyt [J8]. 


Wir zeigen weiterhin 


Lemma 7. Es seien die Voraussetzungen von Lemma © erfiillt und es gelte 


auperdem mit M*(x’, x7, 3) = M(x’, 3) —Me (ev, a) 


(5) pdt (ees 3) =0(, 7 eae | 


real le 
una 
bat adi ex Coat! o (le adie 
ae ee no Ole 
fiir 
(4.14) Ses B= |S Fle 5. 


Dann erfiillt das Feld M gq eine Hélder-Bedingung der Form 


(4.15) [Mg (x’,o) —Mag(e’, »)| = o( | yor log rEral 


fiir 
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Beweis. Die Punkte r’, r’’€ F und y €G seien gegeben und es gelte [e’—2"|S am 
Wir zerlegen die Flache F wie in Lemma 6 in drei Teile A, & und fy, wobei 
die Einteilung beziiglich der Punkte x’ und  getroffen werden soll. Dann gilt 
unter Beriicksichtigung von Gl. (4.9) und Voraussetzung (5) 


1 
= gecoehl see aF. 
ole : 7 |e’—a] “|e” —al -[a—ol* ) 


(4.16) =o yy —e"| fives oars] 


Vina at 


| M*(x’, 2”, 3) 93, 9) a, 
EL 


= o(leoel 


\ ia 
Ferner ist 


| M*(e’, 2”, 3) 9 (3, o) F, 
F, 


(47) =m@'se", 0’) [96.0 dF — f fuer’. 2”, v') mre 2” 196.9) a 
FR, F, 


SAS ls’—r"l £ la—'l —vnlleae4 
=o( BFL Joona, Os Saar? o 7 ). 


Voraussetzung (6) laBt sich hierbei anwenden, da fiir 3c¢/, y’Ch 


lx’ —e"|Sdle’—y'|7 und [g—y'|S3|e’—5 |” 
ist. 
SchlieBlich ist 


| JR a) 860 a =O [Heat brane in) 


x, 


(4.18) 


=ofz Hog — Y L 
Damit ist Lemma 7 bewiesen. 


Als Korollar zu Lemma 7 bemerken wir, da8 im Falle O<a<1 das Feld Mg 
hdlderstetig ist mit einem von y unabhangigen Koeffizienten. Es ist namlich in 
diesem Fall 1- 


jf a2", 8) 969) aF, =O(|x’—2"[® 


wobei diese Abschatzung gleichmiibig fiir alle y €G gilt. 


Wir betrachten nunmehr die in Gl. (4.6) gegebene Kernmatrix und zeigen, 
da8B sie die Voraussetzungen von Lemma 6 und 7 erfiillt. 


Lemma 8. Die Kernmatrizen 
M(x’, 3) = [n(v') VY, Oe, x')] € 
Mes (v’, 3) =V, OG, x’) [(n(v’) x(a) xn) ]~ 


erfiillen beide die Voraussetzungen von Lemma 6 und 7 und zwar bei rein imaginarem 
k=11 gleichmafig fiir alle t=0. 


und 
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Beweis. Es ist 


und 
Ss P62) Ine) xnG)=0( 42). 


Beide Abschatzungen gelten bei rein imagindrem #=it gleichmabBig fiir alle 
t20. Daraus folgt Voraussetzung (3). Nach dem Mittelwertsatz der Differential- 
rechnung ist mit geeignetem # (0<#<1) und mit 3*=34+0(p’ —3) 


(4.24) mle) {ZyP — Vy P} = (y! — a) yr gee O42). 


Hierbei bezeichnet ae 


die Ableitung in Richtung der Normalen im Punkte yx’. 
: 


Setzen wir r=z' —3*, r=|r], v=—r und g(r) = abe: 


, so ist 


n(x’) V2 — Vy ®} 
(4.22) ih Oe 
a= (é 

als dy 


) £05’ — a) x] fn’) 9] + + 5 [(y’ — 3) ne]. 


Fir |y’—3/S3|z’—y|| ist r23|x’—y'| und 


Aas) t de} —o( 4) -0( tr), 
sowie 7 : ; 
(4.24) Ir (5 a} =0(3)= °(paye) 


wobei diese ate ae epee bei rein ana #=17 gleichmabig fiir alle r=0 


gelten. Wegen | n(r’)t9| =O (7) und | n(v’) (y’ —3)|=O(|z’ —p'| - |p’ —3]) erhalten 
wir daher insgesamt 
(4.25) [nV ® — Py} =0( LS), 


d.h. der Kern 3, (x’, 3) erfiillt Voraussetzung (4). 
Um diese Voraussetzung auch fiir den Kern J, nachzuweisen, schreiben wir 


Me (t's 8) — Mee’, y’) = (Vj B — Vy ®) [{n(e’) xn (y’)} xn (1) J~ + 
+ V; O[(n (x) 1()) 2) — (2@) (y’)) ny’), 


und schatzen die beiden Summanden einzeln ab. Es ist fiir |3 —y’|<3}|x’ —’| 


(4.26) 


(4.27) 17, 26,2") — Fy M19’, 2)] =0 (TPF) 
und 
(4.28) [1 (3) (a (z’) Hi) — n(y’) (ne) 2(y’))| = O(| Ae 
so daB wir insgesamt 

' Bat | = 
(4.29) [Me (e’,2) — Male’, 9’)| =O S| 


erhalten. Diese Abschatzung gilt wieder gleichmabig fiir t= 0. 
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Auf entsprechende Weise zeigt man, da die beiden Kerne und Nt, auch 
die Voraussetzung (5) erfiillen, wobei auch hier die GleichmaBigkeitsaussage bei 
rein imagindrem &# gilt. Wir wenden uns nunmehr dem Beweis von Voraus- 
setzung (6) zu. Fithren wir zur Abkiirzung die Funktion 


tev 5p =10)V0G 1) — net PGs) 


ein, so haben wir die Differenz A=f(r’, x”, 3) f(z’, x’, y’) abzuschatzen, und 
zwar unter den Voraussetzungen 


/ 


Selr— hy 


lr’ — 2" |S4le—y'|, |a— 


Mit geeignetem # (0<#<1) und 3*=3+0(y’ —3) erhalten wir nach dem Mittel- 
wertsatz 


(4.30) A= (5 — 9’) Var f(U', 2", 8*)- 
Setzen wir noch 1,=3*—2z’, t,=3*—yx’’ sowie |r,|—7;, —ty=% (j=1, 2) und 
g=—e", so ist d=A,+ A, mit 


(4.31) A, = (5-9) (ey + SE) — Gy) ne) 4 
und : 
(9) eb») 469) 4.2 $F) alowed gL Se) 


Beriicksichtigen wir, daB 


1 dgy(n) 1 dp(72) 


Ys ats Vx ate 


(4.33) 


2* o [em 2” | 
[z’—y =) 
ist, und zwar gleichmaBig bei rein imaginarem 4, und daB ferner 


(4.34) Lacey) p= Ol Taree: ) 


ist, so crhalten wir fiir A, die Abschatzung 


‘|e — 9! 


und zwar gleichmabig fiir t= 0. 
Entsprechend ist wegen 
aa iad od, ad Hiedowd |x’—x’’ | 
A.36 ay fase Sie Olle 8 
Cy, dy, f me) ante ae) Nac 


und 


(4:37) (6 — 9') ft (n@’) 1) — tan") 18)}| = O([2’ — 2] -|8 — 9] 2’ — 0) 
auch 
(4.38) |4,|=0(#oe lero) 
[z’—y’| 
Damit ist gezeigt, daB I, die Voraussetzung (5) erfiillt. Wir haben diese Voraus- 
setzung schlieBlich noch fiir den Kern 9, nachzuweisen. 
Setzen wir zur Abkiirzung fiir die Vektoren 


YyPG Er) =a, VG r")=a,, Vy Py’) =a, Vy P(y',2") =ay, 
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und entsprechend fiir 


[m(x’) xn (g)] Xn (3) =b,,..-, [n(e”) xn(y’)] xn (y’) = by, 
so ist 


(4.39) IMF (u’, v4) — MF (x, 2, o’) =a, WT — ay be” — a3 BF’ + a4 OF. 


Wir zerlegen diesen Ausdruck in die Summanden 


a, (b, — bp — bg + by)~ + (a, — ag — ag + ay) 63° + 
Ga Me) (Pg Vs) g(a; a>) (boy) 


Nach Gl. (4.27) ergibt sich fiir die letzten beiden Summanden sogleich die ge- 
wiinschte Abschatzung. Ferner ist 


|b, — by — bs + ba| =O (|z’— 2” -|3— 


f= VA PAR 
Ja, a, — ag +04, =0(2 Ee). 


) 


und 


Damit ergibt sich insgesamt 
(4.40) | mez (a eG 3) ans me eee ee y’)| =) (oe : 
Da die Abschatzungen bei rein imagindrem #= /t auch gleichmaBig gelten fiir 
alle t=O, ist damit Lemma 8 bewiesen. 

Das Flachenfeld q, (r’, )) =2n(x’) x Vy Or’, y) xa], das wir zur Erzeugung des 
Tensors 8,5 a verwendet hatten, erfiillt die Voraussetzungen in Lemma 6 und 7 
mit «= 2, denn es ist fiir beliebige r’ CF, CG 


(4.41) | Ge (z’, y)| =O beenal 
und } 
(4.42) J V, (3, ») aF| =0(1), 


wobei beide Abschatzungen bei rein imaginarem # gleichmabig gelten. Infolge- 
dessen ist 


1 
(4.43) |M92| =O (eaaP 
(4.44) |Z? g2| =O (log |x’ — y|-*) 
(4.45) |M%g.| =O(1). 


Alle Felder erfiillen nach Lemma 7 eine Hélder-Bedingung. Insbesondere ist das 
Feld M?q, nach dem Korollar zu Lemma 7 hdlderstetig, und zwar gleichmaBig 
fiir alle oe Fiir die Lésung j, der Integralgleichung (4.1) erhalten wir daher 
nach Gl. (4.4) 

(4.46) fo = 92+ Mg, + M?g. 4 fa. 


Dabei ist | f3] =O (1) und erfiillt in r’ eine Hélder-Bedingung der Form O(|r’ —x'’|?), 
gleichmafig fiir alle y. Bei rein imaginarem 4= ¢t konvergiert auch die Neumann- 
sche Reihe der Integralgleichung fiir geniigend groBe t=Ty) (tT) hdngt nur von 
dem Gebiet G ab). Die Abschadtzungen fiir 2 gelten daher in diesem Fall auch 
gleichmabig fiir alle t= 7). 
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Nun wird der Tensor ®,, a durch das Flachenfeld {, (x’, )) gemaB Gl. (3.9) durch 


(4.47) &2(,9)a=— | fa(5, 9) x, OG.2) a 
F 


erzeugt. Wir kénnen Felder dieser Form mittels einer Matrix %t(3, x) als 
J RG. 2) fa 0) ah 


darstellen. In unserem Fall ist 9t(3, r) gegeben durch 
(4.48) N (sz) b= _V, Bl;, x) xb, 


mit beliebigem konstantem Vektor 6. Wir wollen zeigen, daB die Tensoren ®;, 
nur dann singular werden, wenn x und y gegen denselben Randpunkt auf F streben. 
Dazu benutzen wir den Begriff des ,,Lichtweges“ zwischen den Punkten x und 9, 
das ist anschaulich der Weg, den ein Lichtstrahl bei einmaliger Reflektion an F 
zwischen x und y zuriickzulegen hat. Wir definieren diesen Lichtweg, den wir 
mit R(x, y) bezeichnen, durch 


(4.49) R(g, 9) = inf (|x —3| +|9 —a)). 


Wir werden die Tensoren durch geeignete Potenzen des Lichtweges abschatzen. 
Abschatzungen dieser Art sind frither von H. WeEytr [18] gegeben worden. 

Bei dem Beweis des folgenden Lemmas zerlegen wir zu gegebenen Punkten 
r€G, »€G die Oberflache F in drei Teile R, K und F,, und zwar 


F= {3 [¢ -3| S38} 
F, = {,|9 —3| S4$R} 
Fy = §3,5€F, 3¢5}. 


Die Punktmengen / und F, kénnen dabei leer sein. Es gilt nun 


Lemma 9. Das Feld (3, )) und die Matrix N(3, x) seien definiert fiir 3€F; 
t, NEG und stetig fiir 3+: bzw. 3-+r. Es sei ferner 


(1) |9(@5)|=O(|s—5|-*), O<aSa2, 
2) | foG.n)4%|—O(1) (falls «= 2 ist), 


3) |9@.5)—s@’,y)|=0 (eee 
fre Se oy, 

(4) |G, x)| =O(|3 —z|-2) 

(5) JX dF] =0(1) 
und 

(6) [NG x) — Rly’, v)| =o/eet 
fir |g —'|S2|¢—y'. 
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Dann gilt das fiir Feld 
§@, 9) =JRGx) 96, 0) aF, 

die Abschatzung 

[h(z, »)| =O (a 
wober R(x, y) der nach Gl. (4.49) definierte Lichtweg zwischen x und y ist. 

Beweis. Wir setzen 

(4.50) H(t, 9) a NgdF (j=1, 2,3) 

j 
und schatzen zunichst }, ab. Auf F, ist analog zu Gl. (4.9) 


4.54 Ee erly, 
daher 


ISs| = Alles y|* |a— aE ar) 
=9( f aes) - Cr): 
|a—z| R 


Es sei y’ der FuS8punkt des Lotes von y auf F, y’€F falls £, nicht leer ist. 
Fir 3€F, ist |3—r|2$R und fiir zwei beliebige Punkte 3,€h, 3.€F, ist 
la —3e|S2|t—4e|, 1 |3—y'|<4|r—h'|. Unter Verwendung der 
Voraussetzungen (2) und (5) erhalten wir daher fiir 


(4.53) ho = Ny’, v) J (3,0) 4k ae [N(p’, x) — NG, x)] 9») aH, 


Ah. 
, ary 
[Sol = O( )-o( ( la : pth) 
, = tay) 
(4.54) eee ia Lean 
a 1 
=0(7a)- 


Entsprechend erhalten wir fiir }, unter Verwendung der Voraussetzung (3) und (5) 


(4.52) 


y) dk, 


(4.55) h= ( f N (5, z) dF, g(t’, 9) as! N (3, x) [9 (3, b) — ge’, o) aks 


d.h. 
1 s—t'| 
| 93 = 0( 2) +0( lereuafesyae aK) 
ro 


ol) 


Damit ist Lemma 9 bewiesen. 
Als Korollar zu Lemma 9 bemerken wir: Erfiillt q eine Hélder-Bedingung der 
Form 


(4.57) jal, 9)| — 6, 9)|=0( 188) tog BP) 


(4.56) 
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so gilt 
: 1 1 
(4.58) [5(z, »)| =O (ze log +). 


Ist q(x’, y) in bezug auf x’ gleichmaBig hdlderstetig fiir alle yCG, so ist H(z, y) 
stetig als Funktion von y und 4 in G. 

Die in Gl. (4.48) definierte Matrix %(3, x) erfiillt nach Gl. (4.27) und (4.42) 
die Voraussetzungen von Lemma 9. Auch das Flachenfeld g,(z’, y) erftillt nach 
Gl. (4.27), (4.41) und (4.42) die Voraussetzungen von Lemma 9 mit «= 2. SchlieB- 
lich erfiillen die transformierten Flachenfelder Mq,, M?q,, M*g, die Voraus- 
setzungen von Lemma 9 mit jeweils niedrigeren Singularitaten. Die Felder er- 
fiillen eine Hélder-Bedingung nach Lemma 7. Daher ergibt sich insgesamt: 


Satz 6. Es set R(x, 9) der Lichtweg zwischen x und \y, dann tst 


(4.59) [12 (x, 9)| =O ies 


oder genauer 


(4.60) ®9(,)a=— 3" f 92(5,0) x 4G,x) 4 + O( ; log 2). 
FE 


Entsprechend gilt (wegen der Symmetriebextehung ) 
1 
(4.61) | ®21 (, 9)| =O (2). 


Die Abschdtzungen gelten bei rein vmagindren k= ¢t gleichmaPig fiir alle t=%,, 
wobet Tt, nur von G abhdnet. 


Im Fall der Tensoren ®,, und ®,, lassen sich Lemma 6, 7 und 9 nicht ohne 
weiteres anwenden, da z.B. der Anteil 


af (@’, 9) =— =_ ne’) x Hy (a, ®(’, »)) 


des Flachenfeldes 4, (r’, )) von der Ordnung 
hierzu 


x’ —y|~* singular wird. Wir zeigen 
Satz 7. Es ser R=R(v, y) der Lichtweg zwischen xy und y. Dann ist 


1 - 
(4.62) Ry1(t, )) a= — seas af ar (t,o) x V, Oi, x) adF,+0O ee log a : 
F 


Sy es ie d(aK®(, y) 1 1 

(4.63) @eo(z, y) a sean} oe ») V; P(5, x) aF, + 0 (py log =) 
F 

und 


(4.64) |; (¢, 9)| = 0 a (f= 4.9). 

Beweis. Der Tensor ®,, wird durch das Flachenfeld g, in derselben Weise 
erzeugt wie der Tensor St, durch das Flachenfeld q,. Die wiederholte Anwendung 
von Lemma 6, 7 und 9 scheitert aber zuniichst an der Singularitat des Anteils 
qi (x’, »). Wir betrachten daher 


(4.65) M o¥ = f Mr’, 3) 9* Gy) dH, 
its 
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und wollen zeigen, daB dieses Feld die Voraussetzungen (1) und (2*) von Lemma 6 
mit «2 und eine Holder-Bedingung gemaB Lemma 7 erfiillt. Wie in Lemma 6 
und 7 zerlegen wir F in &, # und Fy. Fir die entsprechenden Integrale iiber 
F, und F; erhalten wir (7=|z’ — p]) 


o(-) bzw. 0 (BF hog ar 


da die Abschatzungen dort auch fiir beliebige « > 0 giiltig bleiben. Es verbleibt 
noch das Integral iiber &, d.h. bis auf den konstanten Faktor — a der Ausdruck 


(4.66) (cop) ‘ Ju’) x{[n(a) x V,(aV, 2G, »)) |x OG r)} ab. 


Es ist 
(5 Onl x f 1 (3) XV,(®(, x’) a V, (3, »)) 4 — 
ek = n(e) Fox f (0 ®) (06) x %PG,24) a8, 


Den ersten Anteil formen wir nach dem Satz von STOKES um und erhalten 


f(x’, y) = fe) x (Vy ®(, r’) X ty) (a V, Bg, »)) ds; — 
Fy) 
(4.68) = ap (a V; P) (n’) x uy Plo x) dF, + 


+k f (aV,®) (n@’) xn(3)) OG,x") aR, 
FB, 
Dabei ist © (fj) die Randkurve des Flachensttickes /,, t) der Tangentialvektor 
an ©(F), (t§@=1). Fir 3€F, gilt 
(4.69) [n(e') x (Py BG, 2’) xto)| =O (3). 
Folglich ist 7 


leet) | == o(-5) +0 fap n(x’) lees Vay os r,0)|a¥) 
EO aolseol fetter) 


Ganz entsprechend ergibt sich auch die Giiltigkeit der Hélder-Bedingung auf 
F, und unter Benutzung des Stokesschen Satzes Voraussetzung (2*). Wir kénnen 
daher Lemma 6, 7 und 9 auf Mg} anwenden und erhalten nach den gleichen 
Uberlegungen wie in Satz 6 das Ergebnis 


1 1 
®1(¢, 9) a= — [oF ») XV, PGs, 2) dF + OCF log R): 
F 


SchlieBlich ist der erste Anteil von der Ordnung: O(R~*). Die Anwendung von 
Lemma 9 bietet namlich Schwierigkeiten nur bei dem Integral tiber Fy, die sich 
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aber durch eine entsprechende Anwendung des Stokesschen Satzes umgehen 


lassen. ; 
Wir haben noch den Tensor , zu betrachten. Er wird aus dem Flachen- 


feld j,(x’, ») gemaB ; 
(4.71) Meole,na= tf {oe 6,2) iG) — 5, Woh) BO} ah, 
F 


erzeugt. Der erste Summand des Integranden bereitet keine Schwierigkeiten, er 
ist von der Ordnung O(R7). Bei dem zweiten tritt die Flachendivergenz Vp fy 
auf. Zerlegen wir f, wieder gemaB Gl. (4.2) in 
fo = G2 + Mg, + M?g. + Mg. + fi, 
Vofe=Vode+VoM ge +--+ Vo fi. 
Falls V, q existiert und stetig ist, gilt 

V,(Mg) = f (n(v’) Ve P(e, 2’) Vo, 0) a — Ff (n(e') 9) PAF, 


so ist 


Man erhalt diese Formel, wenn man aus F ein Flachenstiick um x’ ausschneidet, 
von dem Restintegral die Flachendivergenz bildet und nach dem Stokesschen 
Satz umformt. Ferner ist 


(4.73) Vode =— 2-5 (a My ®) + 24° D (ance), 
daher 
(4.74) Voto =Vo do + (M* + M*?4 M*3) Vy go +Vof2+ N** go, 


wobei N** sich aus Produkten von N*, M* und M zusammensetzt. N**q, tragt 
also zu St. nur Glieder von der Ordnung O(R71) bei. {3 geniigt in x’ gleichmaBig 
fiir alle »€G einer Hdlder-Bedingung, dasselbe gilt daher auch fiir die Ablei- 
tungen von f,=MfR. Auf M* Vigo, M*2Vjqg,, M** ogo lassen sich Lemma 6 
und 7 anwenden. Der Kern M*(x’, 3)=n(z’) V,®(3, x’) erfiillt namlich nach 
Lemma 8 die Voraussetzungen von Lemma 6 und 7. Die weiteren Uberlegungen 
verlaufen analog zu denen im Fall des Tensors R®,,, so daB wir auf die Einzel- 
heiten verzichten. 


Wir wenden unsere Ergebnisse an und beweisen: 


Satz 8. In G seien die Vektorfelder j, j' und ihre Divergenzen definiert und 
stetig. Dann erfiillen die Felder 


Ea) = Sri + f Grit Gailey, 
G 


34M 


(4.75) 


ae 


©) = sry gO) + ge f (Grit Geoil ay, 
G 


34 Mb 


mm allen inneren Punkten von G die Maxwellschen Gleichungen 
VX$+dme€ =j 
Vx =o ha= jy 


(4.76) 
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Die Felder €(t), (x) sind im abgeschlossenen Bereich G stetig und es gilt die 
Randbedingung 


(4.77) nx€=o0 auf F. 


Zum Beweis schreiben wir die Felder gema& Lemma1 und 2 in der Form 


C= E+ f Sri t+ Mei) dV = E+ Er 
(4.78) G 
G=GF +3 f Grit Mi) dV = 9" + 9, 
G 
und erhalten nach Lemma 1 zunidchst 


(4.79) VXo* + me =4,— Vxe¥—zon 6* =—j’. 


Da (fiir r€G, y¢F) die kompensierenden Tensoren stetig differenzierbar sind, 
folgt nach Gl. (3.5) ferner 
(4.80) VX O**= vo eG**=7, Vxt* zou S**= 0. 


Aus Gl. (79) und (80) ergibt sich Gl. (76). 


Zum Nachweis der zweiten Behauptung des Satzes zerlegen wir das Gebiet G 
in die Teile G, und G,=G —G,. Zu G, zahlen wir alle Punkte ) €G mit | p 


(c’CF). Wir setzen |y—z'|S = (cr €G) und bilden 


4m G(x) = f Grit Gai) V+ f Git Get) av + 32 


ME 
(4.81) es Go 
SAT, G + 427 & 


Aus der Stetigkeit der Tensoren in G, folgt zunachst 

(4.82) ae n(c ) x6, (z) =. 

Den Anteil ©, schatzen wir ab, wobei wir uns nur mit den wie R-* singular 

werdenden Anteilen zu befassen haben. Nach Lemma 1 und 2 und Satz 6 und 7 ist 
— 82? ¢w en(z’) x &, (zr) 


(4.83) = Jul r') x h(x, 9) d% + 2an(z Sain j) Vj BdF, +O (Je). 


Dabei ist F(G,) aie Randflache des Gebietes wd 

(4.84) he, 0) = J [n@) x Biity) AP) xh Pb, x) ak 
Wir setzen 

H(t, ») = (i(0) %) py) mit 

p(t.) =— f (n@) xV, DG, »)) XV, P(, x) a 


F 


(4.85) 


Es folgt daher 
(4.86) Sh, ») 4V, FE MnO 10 y)) pak, As i) pal. 
Go 
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Nach Lemma 9 ist | p(x, ae d.h. 
(4.87) —822 4 en(t’) XG, (x) =n(e’) x f (nj) (py) + 20%, @) dF, +O (Ve). 


F (Gp) 1 


Dabei verstehen wir unter p(x, )) den Grenzwert, den p(x, y) annimmt, wenn y 


gegen die Randflache F strebt. Wir haben den Ausdruck 
p(t, y) + 2aV,® 


genauer zu diskutieren und kénnen hierbei k=0 setzen. 

Zu y betrachten wir den Spiegelpunkt y*, der durch die Gleichungen 
y=vy’—yn(p’) und y*=y’+yn(y’) (y’CF) fiir geniigend nahe an F gelegene 
Punkte py, d.h. fiir geniigend kleines 9, eindeutig bestimmt ist. Das Feld 

aad 23, : 
p*(¢,y) = 22V, ip" 
erfiillt bei festem ) CG die Differentialgleichungen 


Vep*=0, V,xV,xpr(r,y) =0 


ebenso wie das Feld p(x, y) fir R=0O. Beide Felder sind durch die Tangential- 
komponenten ihrer Randwerte eindeutig bestimmt und nach Lemma 6, 7 und 9 
gilt daher 


, : 
(4.88) p(t, h) = 2a, y*| Bae 


Insgesamt erhalten wir damit aus Gl. (87) 
(4.89) —4aéoen(r’) x E,@)=f (nj) ne’) x(ho gs — fren dF, +0 (Vo). 


: * |e—5*] 
F (Gp) 


Fiir y ¢F verschwindet der Integrand, das Integral erstreckt sich daher lediglich 
tiber den in G gelegenen Teil der Kugelflache |x’ —y| =. Hier ist 


(4.90 "x (V7, Vig= | n0( eee 
ma) 2c Ferrey artaripegy eee 
und wir erhalten die ES Vaan 
(4.94 11 ( —o(le asl 

) |¢@e (eiire= ofeau ; = Ve). 


Ist also |x’ —z| S Je, so ist 
(4.92) | n(x’) x , ()| =O (Ve). 
Aus Gl. (82) und (92) folgt 

lim n(x’) x € (x)= 


ry 
Es ist nach Vorstehendem auch klar, daB das Feld €(r) im abgeschlossenen 
Bereich G stetig ist. Denn bis auf in G stetige Felder ist analog zu Gl. (87) 


EQ) = roe | (i) 22%, 0G, ») +90 0) ae +o 


F 
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Dieses Feld ist aber ebenfalls im abgeschlossenen Bereich G stetig, da fiir y¢€F 


h Gl. (88 
nac (88) | 2a V, P(x, 9) + p(x, »)| = | 


ist. Eine entsprechende Uberlegung zeigt, da® auch das Feld §(r) in G stetig 
ist. Der hieriiber entscheidende Anteil wird durch die Summe der Randintegrale 


SQ) =~ [WK oak, 


42om , 
F 
und A see 
H? (x) = Pa [( n) P(t, y) dF, 
F 


geliefert, wobei 


D(z, v) = ») 7, OG, 2) dE 
F 


ist, und p den Grenzwert fiir )—v’CF bezeichnet. Das erste Integral (5!) 


entsteht nach Lemma 1 und 2 durch Umformung des Tensors 33, wahrend das 
zweite (?) durch eine analoge Umformung wie in Gl. (70) aus dem wie 1/R? 
singular werdenden Anteil des Tensors 8, hervorgeht. Obwohl fiir die Stetigkeit 
von (x) allein die Stetigkeit von j’ nicht ausreicht, ist 


9 | S=s5,, [0 (20%,P +b) dk 


F 


bei stetigem j’ stetig. Wie im Fall des ©-Feldes ist namlich auch hier fiir y’cF: 
~ 1 
| 27 Vy P(x, ) +P (r, y)| = Os) : 
Damit ist Satz 8 bewiesen. 
Fiir die in Satz 5 eingefiihrten ausgearteten Tensoren gilt entsprechend 


Satz 9. In G seien die Vektorfelder j, j' und thre Divergenzen definiert und 
stetig. Setzen wir 1 - 
G=—LiG) + f Gil 9; 0) i(0) a% 
G 


(4.93) 
= gz | Sarl; 0) io) ah, 
G 
1 
C= gy [ S295 90) ah, 
(4.94) G 
ee Ai, 1 ee ore 
pera a ie f Seal 5 0) '(y) Vy, 
G 
so gelten in allen inneren Punkten von G die Gleichungen 
Vx€,—1H=—j, VE,=0 
.96 
4 ? ) Vx De =D. 
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Die Felder ©, und $, (k=1, 2) sind in G stetig und es gelten auf F die Rand- 
eS tx Gi =u x Gyan eS O- 
§ 5. Eigenfrequenzen und Eigenschwingungen 


Eigenschwingungen existieren nach Satz 3 nur, falls w, ¢ und y reell sind. 
Wir kénnen in diesem Fall die Maxwellschen Gleichungen 


(5.1) VxE—twuH=o0, VxH$+iweE=o0 
durch die Substitutionen 

(5.2) S=/eS, C=C, wVfen=o 
zu ; 

(5.3) VxE—4HH=0, Vxh+4HE=o 


normieren, und setzen daher weiterhin e=yw=1 voraus. 

Wir beweisen nunmehr einige allgemeine Eigenschaften itiber Eigenfrequenzen 
und Eigenschwingungen und bemerken zunachst, da8 mit w auch —w ein Eigen- 
wert gleicher Vielfachheit des Randwertproblems 


(5.4) VxE—¢woH=o0, VxH+¢mE€=pd 


nx €=p auf F ist. Denn sind ©, § Eigenschwingungen zum Eigenwert w, so 
sind €, § Eigenschwingungen zum Eigenwert —@. Wir kénnen uns daher in 
der weiteren Diskussion auf positive Werte von w beschranken. Im Fall w=0 
zerfallt das Differentialgleichungssystem. Wir schlieBen diesen Fall zunachst aus. 
Es gilt nun? 
Lemma 10. Exgenschwingungen zu verschiedenen (positiven) Eigenwerten sind 
orthogonal. 


Beweis. w,, @, seien Eigenwerte. Aus den Relationen 


(5.5) Vx, =40,9,, Vx$,=— 40, 
und der Randbedingung n x ©,=0 (y=1, 2) folgt 
(5.6) J V(G,x G,) dV = f (nx &) § dF =0 
G F 
und 
(5.7) [V(G,x,) dV = f ux) §,dF =0, 
G F 
also 
(5.8) o, { ; dV =a, f EG, aV, 
G G 
(5.9) wf, H,.4dV =a, f & Gav. 
G G 


Fiir w} — w3-E0 folgt daraus wie behauptet 


(5.10) {9,9,4dV =[€6,dvV =o. 
G G 


1 Vel. z.B. [5], S. 78 ff. 
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Fiir @,=@, +0 ergibt sich aus Gl. (8), (9) eine Normierungsbeziehung fiir Eigen- 
schwingungen : 


(5.11) [\G|2av =f \gpar. 
G G 


Fir die Eigenschwingungen und Eigenfrequenzen unseres Problems wollen 
wir jetzt eine Integralgleichung aufstellen. Im Fall w=it mit reellem t>=0 
existieren nach Satz 4 und Satz 5 Greensche Tensoren W;,(r, 9; 77). Bilden wir 
mit konstantem a die Felder 


(5.12) Cy = Go(t, 9; Zt) a 
und 
p13) Do = Goi (t, 9; Zt) a, 


so erfiillen sie in G — Kj, die Gleichungen 


(5.14) V.xV_x © + 7? & =o, 
(5.15) V.xV,X 9p +? Bp =o, 


und auf F die Randbedingungen 
(5.16) tx C= 0. Wax )o—D. 


Die Gln. (14) bis (16) bleiben auch fiir t=0 bestehen. Wir nehmen an, sei 
Eigenfrequenz und &, § zugehdérige Eigenschwingungen. Dann gilt 


(5-17) VxE—dwH=d, Vx$+¢oC=od, 


und nx G=p auf &. Wir betrachten den Ausdruck 


(5.18)  f (G-Uxhx€ —E-V,xV,x ©) dV, = (w+ 7) f (G ® aV,. 


of - 
G—Ky G—Ky 


Nach dem GauBschen Satz folgt bei Beachtung der Randbedingungen 


(5.19) nx€=nxG& =o 

aut f ; 

(5.20) (w2 + 72) f (€, G) dV, = f n[ExV,xG—& xV,x €] dF,. 
G-Ky Fi 


Lassen wir jetzt 7—>0 streben und beriicksichtigen nach Lemma 3 und 4 die 
Singularitaten des Feldes © 


(5.21) C= ox 

so erhalten wir 

(5.22) (w? + 7) f (Gp ©) dV, =— 4a 2wa H(b). 
G 


In Gl. (22) ist a ein beliebiger konstanter Vektor, wir kénnen sie daher auch 
in der Form 


(5.23) (wo? + asl Gia (z, 9; 47) E(t) dV, =—4a¢oH(y), 
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bzw. unter Benutzung von Satz 1 oder Satz 5 in der Form 

(5.24) (w? + 7) f Gai (0,25 47) EQ) dV, = 402 Hb) 
G 

schreiben. Analog ergibt sich aus dem Ansatz 


(5.25) f (Go %xUxS— 9 -VxXV,x Go) dV; = (o? +7) f (Go 8) a, 


G—Ky G—Ky 
bei Beriicksichtigung der Randbedingungen 


nxXV,XxH$=nxV,xXHo>=o0 auf F 


und der Singularitaiten des Feldes 


yo =V,Oxa+-- 
die Relation 
(5.26) (0? + 7?) f Gre (y, £5 47) H(z) dV, = 4 Zw E(p). 
G 
Setzen wir 
ea ae © — G,2(t, 9; a8) 
(5.27) U(r, 95 ¢t) = tle. ae 5 


wobei © die dreireihige Nullmatrix ist, und fassen die Eigenschwingungen gemaB 


EG) 
.28 = 
(5.28) PO =( 
zu einem sechsdimensionalen Vektor zusammen, so kénnen wir fiir Gl. (24) und (26) 
2 2 
(5.29) OEE [ 86.9: 47) (0) ah =PO) 
G 


schreiben. Jede Eigenschwingung des Randwertproblems — nach Gl. (28) zu- 
sammengefaBt — ist eine Lésung dieser Integralgleichung. Nach Satz 1, 2 und 5 
ist der Kern &(r, t); ¢t) reell und symmetrisch 


(5.30) QU" peer) (ee 
seine Elemente sind von der Ordnung oar}: Der iterierte Kern ist daher 


ebenfalls reell, symmetrisch, quadratintegrabel und nicht ausgeartet. Nach einer 
bekannten SchluBweise von HILBERT [8], S. 69, besitzt daher die Integralgleichung 
(29) abzahlbar unendlich viele Eigenwerte und zugehérige Eigenlésungen. Wir 
haben noch zu zeigen, daB jede Eigenlésung der Integralgleichung (29) auch 
Eigenschwingung unseres Randwertproblems ist. 


Wir beweisen zunidchst 


Lemina 11. Der Kern 
(5234) Yo(t, 0; ¢t) = f Be, 4; 0) 2G, 9; 42) AV, 
G 


ist reell, symmetrisch und quadratintegrabel, 
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Da &(z, 3; 0) und (3, ); ¢7) reell und von der Ordnung 
of | eer | (See 
nal) ay (=F) 


sind, ist der Kern &, reell und von der Ordnung o(-+, d.h. quadrat- 
integrabel. Aus Gl. (30) folgt weiterhin lz—9 


(5.32) Le (t, 9; 47) =J& (v.43 ¢7) &(3, 4; 0) dV, 

so daB wir zum Nachweis der Symmetrie lediglich die Relation 
(5.33) [&os ét) 8,9; 0)aY=J&ea: 0) 2, 9; 47) dV, 
bzw. 


(5.34) Jub s: éT) Goi (5, 93 0) ee! G12 (t, 3; 0) Gai (3, 9; ¢t) AV, 
und 
(5.35) J Sales: £T) ®o(3, 93 0) jog Goi (t, 33 0) Grog, 9; 7) AV, 


nachzuweisen haben. 


Es sei f(r) ein in G stetig differenzierbares Feld. Setzen wir zur Abkiirzung 
(20) 


(5.36) Fe f Bin 0; ¢2) FW) 4% = Oh, 
und bilden die Felder . 
©. = Gif, OF = O81} 
Ga-Lf+Ohf, G=-Li+Ghi 
e°” E, = Ge Ha, GF = Oi. SF 
S%i=—7%2.+G%b, Sf=— ) OF + Ge OF, 
so ist Gl. (34) gleichbedeutend mit ©,= G7. Nach Satz 8 und 9 gilt aber in G 
(5.38) Vx$.—1 & =f, Vx or + 26h =I, 
(5.39) Vx G+ H.=o, Vx GE =; 
(5.40) Vx6,—¢$,=—%, Vx6t+r$*=— GF, 
(5.41) Vx D1 =—05 Vx $i —r Gf =p, 


sowie VG,=V Hz=0. Auf F gelten die Randbedingungen 
t~Ganx Cc =p (=472) und 
We — 1 oy 0: 


Nach Gl. (39) und (40) ist V$,=V,=0, also ,= 0. Eliminieren wir nun aus 
Gl. (37 bis 41) die ttbrigen Felder, so erhalten wir fiir ©, und Gf 


V5xG + ?VxX& =—Vxi 
(5.43) Vax GF + 2Vx GF =—Vxf 


(5.42) 
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sowie! 
(5.44) nxV,x(G,— Gf)=0 auf F. 
Aus Gl. (42 bis 44) folgt unmittelbar ©, = G* und damit die Behauptung Gl. (34). 


Ganz entsprechend beweist man Gl. (35). Wir haben bis jetzt gezeigt, daB 
die Eigenschwingungen Lésungen der Integralgleichung 


(5.45) ott) 5 (r,.9; 47) p(y) dV, = (2) 


sind, deren Kern reell, symmetrisch und quadratintegrabel ist. Wir zeigen nun 


Lemma 12. Jede Eigenlisung der Integralgleichung (45) lefert eine Evgen- 
schwingung des Randwertproblems (5.4). 
s 
Beweis. Es ser p= (*) Eigenlésung, d.h. es gelte — mit normiertem © 
und § — 1 
—o 8, 63,E=6 


(5.46) 
—= 6 S$, OFE oe = §. 


Dann gibt es auch Zahlen o, und o, mit o, o,.=«a, so daB 


o, 8, E=1H 


(5.47) 
— 405 OMe D = 6 


ist ([8], S. 69). 
Aus Gl. (47) lesen wir ab, daB die Divergenz der stetigen Felder &, § in G 
verschwindet?. Bilden wir daher 
G, SS £09 Po s 
De = — ¢0y(— 30 + B22 9) 
(5.48) 


= 0 ( Bye Gis 6) 
OT 
D, = — ¢o, 63, &, 
so ist C=, H=H, und nx G& =n x €,=0 auf F. Ferner gilt nach Satz 8 und 9 


VxG@—¢ D275, Coun VE, =0 


Vator = 
(5.49) e : 
Vx €,+7$,=0 
VX 9, — t © =— 40, G,. 


Wegen V$,=V,=0 und n$,=0 auf F folgt zunachst wieder §,=o, dh. 
Vx &, i C09 Dy = 0 

(5.50) ~ te 0h ae) 
Vx (Vx, L = &, =D. 


- - 2 


1V, x E (n=1, 2,...) soll die n-malige Anwendung d 
5D, a g des Operators rot andeuten. 
2° Vel. [17], Lemma 5. ‘ | 
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Da ferner in G V €,=0 und auf F auch nxVx,=0 ist, erhalten wir schlieBlich 


Vx ¢,— 409 Soy =D 
(5.51) 


a, re 
Boo Civ; <= 0. aul PF. 
O2 


Bei nicht identisch verschwindenden Feldern €,, , kann dieses Gleichungssystem 
aber nur bestehen, wenn o —t?=q? fiir eine Eigenfrequenz w ist. 

Wegen der Normierung ist ferner o,=o0=-+-|o —7?, d.h. &, §, sind Eigen- 
schwingungen zur Eigenfrequenz m. Da der — positiv definite — Kern &, (r, ); g 7) 
quadratintegrabel ist, gilt schlieBlich ein Entwicklungssatz fiir die durch ihn 
quellenmaBig dargestellten Felder. Wir formulieren unsere Ergebnisse in folgen- 
dem 


Satz 10. Den Eigenlésungen der Integralgleichung 
(wo? + za L(t, 9; 47) p(y) dy =p) 


mut reell-symmetrischem, quadratintegrablem und positiv definitem Kern entsprechen 
umkehrbar eindeutig die Eigenschwingungen des Randwertproblems. Es gibt eine 
abzahlbar unendliche Folge von Eigenfrequenzen w,,==0, die sich im Endlichen nicht 
hiuft. 

Jeder Eigenfrequenz w,, entspricht eine Eigenfrequenz —o,, gleicher Vielfachheit, 
wobet die Eigenschwingungen zur Eigenfrequenz —w,, konjugiert komplex zu denen 
von w,, sind. 

Jedes Feld § (x), das mittels des quadratintegrablen Feldes q durch 


f(x) eu @ (x, y; #7) q(y) dV, 


quellenmapig darstellbar ist, kann in eine gleichmadfig konvergente Rethe nach den 
Exgenlosungen ,, entwickelt werden. 


§ 6. Asymptotische Gesetze fiir Eigenschwingungen und Eigenfrequenzen 
Wir wollen in diesem Paragraphen asymptotische Gesetze fiir die Eigenschwin- 
gungen und Eigenfrequenzen der Hohlraumstrahlung herleiten. Dazu unter- 
suchen wir das Verhalten der Tensoren ,4(r, ; 77) und G., (x, ); 77) fiir r>—&, 
aus dem wir auf asymptotische Gesetze fiir die Eigenschwingungen und Eigen- 
frequenzen schlieBen kénnen. Setzen wir zur Abkiirzung 


(6.1) Gre, 9; 4t) =— ee | Gy 2(, 3; 0) G13, 95 47) dV, 
G 
(6.2) Gye, 9; ¢t) = — Gar | Sunt 3; 0) G23, 9; 47) dV, 
so ist G 
G, O 
(6.3) 82(6,95 49) =( < i 


und aus Satz 9 folgen in bekannter Weise! Entwicklungssadtze fiir die Tensoren 
Us Or und Oy. 
1 Vgl. [3], S. 114ff., [16], S. 227ff. 
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Zur Vereinfachung der Schreibweise wollen wir dabei wieder 


_ { ,(e) 
pn (t) = Le a 


setzen, also unter ),,(r) einen 6-dimensionalen Vektor verstehen, der sich aus 
den Eigenschwingungen G,,(r), ,(z) in der angegebenen Weise zusammensetzt. 
Ferner verstehen wir unter p,,(r) )> (y) die 6-reihige Matrix, die sich ergibt, wenn 
man die Spalte ,,(r) mit der Zeile p> (y) nach den Regeln der Matrizenrechnung 
multipliziert. Dann gilt 


Satz 11. Sind &,,(r), 5, (x) die orthonormerten Eigenschwingungen zur Frequenz 


@, ==0, und st 
G,, 
m=( eg): 


so gilt 
2 (oe) 
(6.4) (e950) - Seo. = > yi See eee 
oder — damut gleichbedeutend — 
Ritne te sv Gn Ene) Gi (vy, Si ( )) 
(655) Gels, 9; 0) — Gee, 9; 7) = 71D’ Re 
(6.6) Gy (t, 9; 0) — Gye ys 47) = 7? ae 2210 95 (0) 


wobet > bedeutet, dak nur tiber die positiven Eigenfrequenzen w,, zu summieren 
ast. Die Rethen Gl. (4), (5), (6) konvergieren fiir alle xy und y aus G gleichmafig. 

Wir wollen nunmehr auf beiden Seiten der Gl. (5) und (6) die Spur (also die 
Summe der Hauptdiagonalelemente) bilden und speziell das Verhalten der Spur 
fiir t—>co untersuchen. Wir zeigen zundchst 


Lemma 13. Es ist 


, : ; ise De Rg alee 1 
(6.7) Hag SPU (Ge (zr, ); 0) — Getz, 5 4%)) ae | Oar 

(6 : kane. Tt 1 
(6.8) a Spur (Gy (t,); 0) — Ga (x, 9; ¢7)) = oe silbe4h 


wober R, die kiirzeste Entfernung von y zu F ist und die Abschatzung gleichmaBig 
Henge alley T= To gilt. 


Beweis. Beriicksichtigen wir nach Gl. (2.3) und (3.32) die Singularitaten der 
Tensoren ©,, und ®,,, so erhalten wir 


Spur { ®12(6 350) O16, 91¢7) ah =—2 [4 7,06, 0) ar,4 


(COs $ 
O See be A Mi 21 D ees ,¢ dave ee 
=i, / Pi2(t, 3; 0)|- [Gor 9 nah) +f fib )| [g—t]? Mm), 


1g—z| 
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Den ersten Summanden formen wir nach dem Gaufschen Integralsatz um und 
erhalten 


é 1 
iar faa] P60) = 4x Oley + [ ee 
(6.10) = 47 Gz, ») oC ieee os ma an) 

1 
= 4x ®(r, ») +0 (Zz). 
Nach Satz 6 gilt gleichmabig fiir alle r=t,) 
~ 1 
(6.11) |B.166.95 <2)| =0(ger), 
wobei R,, der Lichtweg zwischen 3 und y ist. Da R,,=|3—y| ist, gilt ferner 
gleichmaBig fiir alle t= T, 
1 
(6.12) |Go1 (95 ¢7)| =O (react 


Die weiteren Summanden in Gl. (9) kénnen also durch 


dV, AV; 
of. ne o( aca) 


abgeschatzt werden. Da R,,=R, und R,y2R, ist, kénnen diese Ausdriicke 


auch durch 
o( 1 {i av, al 1 
Rie lecelhlesol? Ue 


O f ays )=0 : 
ce la—o|? |a—z]? (aH) 


ersetzt werden. Insgesamt ergibt sich nach Gl. (1) und (9) 


Spur (Gz (z, 9; 0) — Ge (x, 9; ¢7)) 


(6.13) = el rer == D(,, »)] Oly i st 


Fiihren wir jetzt den Grenziibergang ) x aus, so folgt 


(6.14) lim Spur (Gee, 9; 0) — Ge, 9; ¢x)) =F + O(S4), 


wobei die Abschatzung O (1 /R}) gleichmaBig beziiglich aller t= Ty gilt. Die zweite 
Behauptung von Lemma 13 1aBt sich wortlich ebenso beweisen. 


Da aun Spur G, (x) GF (9) = GF @) E, (0) 


a Spur 6, (2) Ee) = |G, Ol? 
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ist, erhalten wir aus Gl. (5) und (6) folgende asymptotische Relationen fiir T—> co 
und r€G, r¢F: 


Ile 1 

(6.15) ds, “we (2 + 7) Pyne & 
| n(x)? mer 

(6.16) Der re T?) Ont * 


Wir wollen zeigen, daB diese asymptotischen Relationen tiber G integriert werden 
diirfen: 
Lemma 14. 
1 V 
(6.17) > w2 (w2 + 1?) Wy 


On > 


wobet V das Volumen von G ist. 


Beweis. Es sei G, der Teilbereich von G, fiir dessen Punkte r€G5: R,2 0 ist. 
Integrieren wir die Relation 


En (x 1 i) 1! 
(6.18) DH =. ake salle 


zunichst tiber G,, so folgt wegen der gleichma8igen Konvergenz der Reihe fiir 
alle r€G und alle t=0 


E , 1 V 
an wo? (wp + T?) 27 
(6.19) ©, 16, ()|2 ; : ve 
= 4 <M \E)\ palontas men Va 
ail > are 4 fa ++ foyer. 
G-G5"— G—G5 Gs z 
Nun ist ; 
fav =0(0) und f( ;)4%=0(4) 
G—Gy5 G5 Ry/ 


d.h. insgesamt 


(6.20) 


fee) 
, V 
os ay 
On ress 20 


n=1 


Setzen wir daher t= ay SOMSE TU - ONG ——oo 
bt 


(6.21) 


co 

«> 
— we (wz e 1?) 22 
was zu bewelsen war. 


Wie schon T. CARLEMAN [71], [2] zeigte, folgen nach bekannten Tauberschen 
Sdtzen nunmehr unmittelbar weitere asymptotische Relationen iiber Eigen- 
funktionen und Eigenfrequenzen. Wir benutzen ([2], [7], vgl. auch [9], [14]): 

Lemma 15. Ist a,>0O und 7 oo konvergent, wobei i, >O eine monoton 


n=1°” 


gegen co strebende Folge ist, und gilt mit konstanten H +0 und 0<oe<1 


[o.2) 
(6.22 ee, ee 
) 2 pra Hs fir x—>© 


so folgt 
(6.23) 


wh 
an a(1—@) 
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Wenden wir diesen Tauberschen Satz auf Gl. (15, 16, 17) an, so erhalten wir 


(6.24) 
(6.25) 


(6.26) 


~ Gn) t 
wz ae 


0<an<t 
[Sn(x)[2 ot 
ieee ae 
wz I 
0<on<t 
1 Vit 
Fay rom 
(62) TU 
0<an<t de 


und weiterhin durch partielle Integration 


(6.27) 
(6.28) 


(6.29) 


Gl. (29) ist von H. Weyt [17], [18] auf anderem Wege hergeleitet worden. 


G, (2) 22 ~ 


25 


>») 
a 
0<an<t 3 
2 “8 
» |, (x)| Wineat , 
0<an<t 
7- 
{3 
(oat ald 
3.77 
0<an<l 


Be- 


riicksichtigen wir bei konstanten ¢, w die Substitutionen Gl. (5.2), so kénnen wir 
unsere Ergebnisse in folgendem Satz zusammenfassen: 


Satz 12. Es seven ©, 9, Eigenschwingungen zur Eigenfrequenz w, der Maxwell- 


schen Gleichungen 


VES —Zorh=—o0, VXO+eoeC=o0 mG 


nx €=o auf F mit positiven Konstanten ¢ und uw, die durch 


(6.30) 


+{& 6,,dV =+ [8 Dm dV = Onn 
le — 
G 


orthonormiert sind. Dann gelten folgende asymptotischen Relationen (x CG, tov) 


(6.31) 
(6.32) 
(6.33) 
(6.34) 


(6.35) 


§ a 
eee 


0<an<t On 
Bel x)? wwe ley)s 
On IU 
0<an<t 
ws By C7 
= |G, (t)[2>~ ue (e u)® 3 02 , 
0<yn<t 
8 3 
a ID, (x)|2?~e(ew)i— 2°? 
i a 
0<an<t 
3 
~ je. 7 ls 
1 (emi Ve 


0<an<t 


348 


[79] 


Mu.ier und NreMEYER: Asymptotische Gesetze der Hohlraumschwingungen 


Literatur 


CaRLEMAN, T.: Propriétés asymptotiques des fonctions fondamentales des mem- 
branes vibrantes, C.R. du 8° Congrés des Mathématiciens Scandinaves, Sund 


1935, S. 34—44. 


] CARLEMAN, T.: Uber die Verteilung der Eigenwerte partieller Differentialglei- 


chungen. Ber. Sachs. Akad. d. Wiss. 86, 119 (1936). 


) Courant, R., u. D. Hireert: Methoden der mathematischen Physik, Bd. I, 


Berlin 1931. 


] Courant, R.: Uber die Schwingungen eingespannter Platten. Math. Z. 15, 195 


1922). 

Oe G.: Elektromagnetische Wellenleiter und Hohlraume. Stuttgart 1955. 

Gunter, N. M.: Die Potentialtheorie und ihre Anwendung auf Grundaufgaben 
der mathematischen Physik. Leipzig 1957. 

Harpy, G.H., & J. E. LirttEwoop: Notes on the theory of series (XI): On 
Tauberian theorems. Proc. Lond. Math. Soc., Ser. II 30, 23—37 (1930). 

Hivtpert, D.: Grundziige einer allgemeinen Theorie der linearen Integralglei- 
chungen. Leipzig u. Berlin 1912. 

KarAMATA, J.: Neuer Beweis und Verallgemeinerung der Tauberschen Satze, 
welche die Laplacesche und Stieltjesche Transformation betreffen. J. reine 
angew. Math. 164, 27—39 (1931). 

MULLER, CL.: Zur mathematischen Theorie elektromagnetischer Schwingungen. 
Abh. Dtsch. Akademie der Wiss. zu Berlin, 1945/46, Nr. 3. 


] Mutter, Cx.: Grundprobleme der mathematischen Theorie elektromagnetischer 


Schwingungen. Berlin 1957. 


] MULLER, CL.: Randwertprobleme der Theorie elektromagnetischer Schwingungen. 


Math. Z. 56, 261 ff. (1952). 


] MUxier, Cxi.: Grundoperationen der Vektoranalysis. Math. Ann. 124, 427ff. 


(1952). 


] Prere1, A.: On the eigenvalues and eigenfunctions of elastic plates. Comm. on 


Pure and Appl. Math. 3, 1ff. (1950). 


] PLerjEL, A.: Green’s functions and asymptotic distributions of eigenvalues and 


eigenfunctions. Proc. Symposium spectral theory and differential problems, 
1955, S. 439—454. 
Riesz, F., u. B. Sz. Nacy: Vorlesungen iiber Funktionalanalysis. Berlin 1956. 


] Weyt, H.: Das asymptotische Verteilungsgesetz der Eigenwerte linearer par- 


tieller Differentialgleichungen (mit einer Anwendung auf die Theorie der 
Hohlraumstrahlung). Math. Ann. 71, 441—479 (1912). 


] Weyt, H.: Das asymptotische Verteilungsgesetz der Eigenschwingungen eines 


beliebig gestalteten elastischen Kérpers. Rend. Circolo Mat. Palermo, 39, 
14—50 (19145). 

Wevt, H.: Ramifications, old and new, of the eigenvalue problem. Bull. Amer. 
Math. Soc. 56, 115—139 (1950). 


Institut fiir Reine und Angewandte Mathematik 
Rhein.-Westf. Technische Hochschule 
Aachen 


(Eingegangen am 16. Dezember 1960) 


Fin Existenzsatz fir die Lisungen gewisser Gleichungen 
mit Nebenbedingungen bei beschrankter Nichtlinearitat 


Hans EHRMANN 


Vorgelegt von H. COLLATZ 


Eine Reihe von Existenzsatzen fiir die Lésungen von allgemeinen Randwert- 
aufgaben bei gewissen nichtlinearen Gleichungen geben als wesentliche Bedingung 
fiir die Lésbarkeit die Existenz einer Lésung von einer zugeh6rigen linearen Rand- 
wertaufgabe an, z.B. in Form des Nichtverschwindens einer Determinante. Dies 
ist gewohnlich dann der Fall, wenn die in der Gleichung auftretenden nichtlinearen 
Ausdriicke durch eine von der Losung unabhangige Funktion beschrankt werden 
kénnen. Hierzu gehéren vor allem einige grundlegende Arbeiten von ROBERTO 
Conti! und WiLi1AM M. WuyBurN? tiber die Existenz der Lésung von nicht- 
linearen Differentialgleichungssystemen mit allgemeinen linearen Randbedin- 
gungen. 

Wesentlich in den Beweismethoden ist hierbei die Eigenschaft dieser Glei- 
chungen, daB ihre Lésungen durch endlich viele Zahlenwerte c; (¢=1,..., 2) 
(i.allg. die Anfangswerte) festgelegt werden kénnen, von denen sie stetig abhangen. 
Die allgemeine Lésung l4Bt sich dann in Form einer Gleichung y=F (y(x), ¢;) 
formell angeben. In die Randbedingungen eingefiihrt liefert dies ein Gleichungs- 
system fiir die Zahlenwerte c;, dessen Lésung unter Benutzung der Beschranktheit 
der auftretenden nichtlinearen Ausdriicke mit Hilfe des Brouwerschen Fixpunkt- 
satzes gesichert werden kann. 

Hier wird ein entsprechender Satz fiir gewisse nichtlineare Gleichungen in 
Banach-Raumen bewiesen, der einige dieser Existenzsatze in verschiedenen 
Arbeiten als Spezialfalle enthalt? und die Beweise vereinfacht. Es gelingt da- 
durch, geschlossene Ausdriicke fiir die Bedingungen der Existenz anzugeben, 
sowie gleichzeitig auch mégliche allgemeinere Nebenbedingungen, z.B. gewisse 
nichtlineare Randbedingungen bei nichtlinearen Differentialgleichungssystemen 
mit zu erfassen. 

Die Beschrankung der Nichtlinearitat durch ein festes Element, die in dem 
Satz als wesentliche Voraussetzung auftritt, ist allerdings eine in vielen Fallen 
zu starke Einschrankung, jedoch kann der Satz als Grundlage weiterer Existenz- 


1 R. Conti [4], [5], [6]. 

2 W. M. Wuyzurn [2], [3], s. auch die dort angegebenen weiteren Literaturstellen. 

3 Nicht enthalten sind dagegen die Ergebnisse der Untersuchungen von R. ConrTlI 
[6] iiber gewisse pseudolineare Differentialgleichungssysteme. 
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siitze dienen, wie dies z.B. H. EpHESER [20] in einem entsprechenden Fall bei 
speziellen Randwertaufgaben mit Differentialgleichungen zweiter Ordnung durch- 
fiihrte. 

In einigen Beispielen mit speziellen Gleichungen der betrachteten Art wie 
nichtlinearen Gleichungssystemen, Differenzengleichungen und vor allem Rand- 
wertaufgaben bei gewohnlichen Differentialgleichungen soll die einfache Anwen- 
dung des Satzes gezeigt werden. Insbesondere wird der Nachweis periodischer 
Lésungen der auf hdhere Ordnungen verallgemeinerten Duffingschen Differential- 
gleichung gefihrt. 


1. Wir betrachten Gleichungen der Gestalt 
(1.1) Lu=Mu 


mit einem stetigen linearen Operator L und einem stetigen Operator M. Beide 
Operatoren sind in einem Banach-Raum B, definiert, und ihr Wertebereich liegt 
in einem Banach-Raum b,. 

Die Operatoren L und M besitzen folgende Eigenschaften: 

a) Der Nullraum N von L ist endlich dimensional, d.h. ede Lésung der Gleichung 


Lu=O,  (Nullelement in B,) 
laBt sich in der Gestalt 


n 
U= > by Dy 
a | 


mit n linear unabhangigen Elementen p,€ By, v=1,...,, der Basis von N, dar- 
stellen. 


b) Es extstiert ein vollstetiger linearer Operator L> derart, daB die Gleichung 


Lu=w, wE€B, 


durch 
(4:2) u=Lw 
gelost wird! und die Gleichung 
(1.3) u=>c,g,+L*Mu 
v=1 
fiir jedes feste System der c,, v=1,..., n, eine eindeutige Loésung besitzt. 


c) Ferner existiere eine Zahl m mit 
(1.4) |Mul<m fiir alle wCB,. 
Unter diesen Voraussetzungen gilt der 


Existenzsatz. Es seien F,, w=1,..., 0, stetige lineare Funktionale aus Be 
(dem konjugierten Raum von B,) und H wu» 6=1,...,, Stetige, nicht notwendig 
lineare Funktionale tiber By, mit der Eigenschaft 


(1.5) |1,,[«]|<h — fiir alle w€ By und einer festen Zahl h. 


c ct braucht nicht zu existieren. Die Forderung besagt lediglich, daB eine spezielle 
Losung existiert. 
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Dann besitzt die Aufgabe (1.1) mindestens eine Lésung, die den Nebenbedingungen 


(1.6) El ial, AG cca 10s 
gentigt, wenn die Determinante 

(1.7) . det (F,,[y,]) == 0 

ast. 


2. Zum Beweis treffen wir folgende Verabredung 


Metts (0,(C) 02,0, (c)) bezeichnen wir als Sammelbezeichnung eine vom 
Vektor c=(cy,...,¢,) und eventuell noch anderen Variablen abhingende Vektor- 
funktion mut folgenden Eigenschaften: 


a) Die Komponenten o; hingen stetig von den Komponenten c;, J=1,..., 0, ab. 
B) o «st fiir alle festen Werte der eventuell noch auftretenden anderen V ariablen 
beschrdnkt, d.h. es existiert eine Konstante K = K(a) mit 
ioe CA pny Mod GUE Ce 
Wollen wir hervorheben, dap o nur vom Vektor c abhingt, so schreiben wir o=o%*. 


Man beachte, daB o von Fall zu Fall verschiedene Ausdriicke darstellen kann. 
Zum Beispiel gilt! ¢ -o=o. 

Mit dieser Verabredung V besitzt die Vektorgleichung c=o* nach dem 
Brouwerschen Fixpunktsatz eine Lésung c= (c,), denn der beschrankte n-dimen- 
sionale Bereich |c;| << AK wird durch o* stetig in sich abgebildet?. 


3. Bewets des Existenzsatzes. Nach 1a) und 1b) ist eine Lésung w von (1.3) 
auch Lésung von (1.1). Wir zeigen, daB mit der Verabredung V fiir eine solche 
Lésung 


(3.1) L2Mu=o 
gilt. 

Ware dies nicht der Fall, so gabe es eine Folge c¢y)= (¢;(,)), R=1, 2,..., die 
gegen einen Vektor c konvergiert, wobei gleichzeitig die nach 1b) eindeutig be- 
stimmten Lésungen wu, der Gleichungen 


n ~~ 
Up =>) Coa) Py + LIM uy, pay ne 
ia 


nicht gegen die Lésung w von (1.3) mit c=(c,) konvergieren. Es existierte dann 
eine positive Zahl 6 und eine Teilfolge w,, mit 
(3.2) lle — uw, =6, 7f=1,2,.... 

Wegen der Vollstetigkeit von L> und (1.4) kénnen wir annehmen, daf fiir 


diese Teilfolge bereits LM u,, gegen ein Element @¢€ B, konvergiert. Dann 


konvergiert “,, gegen eine Losung i=) Cy P,+W von (1.3); denn es gelten die 
Ungleichungen ors 


lim fry, — a] Sim ¥ Josey — 66] “Hd + im |Z 2M my, — | = 0 
j— oo = OO 7=1 — oo 


1 Die Funktion o ist also 4hnlich zu handhaben wie das Landausche Symbol O 
in der Zahlen- bzw. Funktionentheorie. 
2 Siehe etwa G. SANSONE e R. Cont! [7], S. 441 oder L. S. PonTRJAGIN [J/], S. 79. 


352 Hans EHRMANN: 


und 


SSM ek Re Renee 


n ~ 
tas 
Up; — DX Co1%) Py L M u;, 
v= 


Dies steht wegen der Eindeutigkeit der Lésung von (1.3) im Widerspruch 
zu (3.2). 

Beachtet man, da® ein stetiges lineares Funktional auch beschrankt ist?, 
so ergibt sich mit (1.5) und (3.1) durch Einsetzen von (1.3) in (1.6) wegen der 
Stetigkeit der Funktionale F, und H,, mit V sofort 


% 6 File] = 0%. 


Wegen des Nichtverschwindens der Determinante (1.7) erhalt man schlieB- 
lich durch Auflésung dieses Gleichungssystems 


C=O, 
woraus nach den obigen Bemerkungen die Existenz einer Lésung folgt. 
Zusatz. Die Existenz einer Lésung von (1.1), (1.6) ast unter den tibrigen Vor- 


aussetzungen offenbar auch dann schon gesichert, wenn wir in 1b) voraussetzen, 
daB es eine Lésungsschar u=u(c) von (1.3) gibt, die stetig vom Vektor c abhanet. 


Vollstetigkeit des Operators L sowie Eindeutigkeit der Losung von (1.3) 
braucht dann nicht gefordert zu werden. Wir haben nur jeweils Lésungen der 
stetigen Schar zu nehmen, fiir die (3.1) gilt. 


Ist die Bedingung (1.7) nicht erfiillt, so hat schon die Aufgabe 
Lu=0,, Eu]=A,, 


Leen iis 


i 


im allgemeinen keine Loésung. Man kann den Satz auch so aussprechen, daB 
unter den tibrigen Voraussetzungen [auBer (1.7)] dann eine Lésung von (1.1), 
(1.6) existiert, wenn die vollhomogene Aufgabe 


Lu=0,, Fluj=0 


nur die Null6sung w=@, besitzt. 


4. Zur Illustration des Satzes mége zunachst der Fall der gewéhnlichen nicht- 
linearen Gleichungssysteme behandelt werden: 


(4.4) Lu=Mu=g(u) 


sel ein nichtlineares Gleichungssystem, L eine mXn-Matrix vom Range m<n, 
u ein n-Vektor (m x 1-Matrix) und g (wv) ein m-Vektor mit der Eigenschaft |g; (w)| Sa, 
fiir alle uv, 1=1,...,m. Ist m<n, so sind die Banach-Raume?2 B, und B, ver- 
schieden (namlich n- und m-Vektoren). Die Bedingungen 1a—c kénnen erfiillt 
werden. Insbesondere existiert wegen des Ranges m von L eine n Xm-Matrix L-. 


1|F[u]|<N|ull, s. z.B. Kormocorov & Fomin [10], S. 77. 
* Man wahle irgendeine der iiblichen Normen, z. B. 


Cenc, n 
lel =[/ Sb oder Di |u;| u.a. 
t=1 t=1 
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Man erhalt diese z.B., indem man 2 —m passende Spalten der Matrix L streicht, 
die gleichzahligen Zeilen von jee gleich Null und die tibrigen Elemente von ibe 
gleich den Elementen der Inversen der reduzierten Matrix L setzt!. Ein Funda- 
mentalsystem y,, v=1,...,2%—m, bekommt man durch Auflosung der Glei- 
chungssysteme Lw—0 mit m—m linear unabhangigen Wertesystemen der freien 
Unbekannten?. 

Es kénnen dann »—m zusatzliche Bedingungen® der Form (1.5), (1.6) vor- 


gegeben werden. Uber die Lésbarkeit des gesamten Systems entscheidet die 
Determinante (1.7). 


5. Einen weiteren Gleichungstyp, auf den unser Satz anwendbar ist, stellen 
die nichtlinearen Differenzengleichung dar, bei denen die Nichtlinearitat in einem 
additiven beschrankten Glied vorkommt. Der Operator L hat in diesem Falle 
die Gestalt 


(5.1) LiE=A" f(x) +a, A" 1f (x) +---+a4,f (x), 20, ganz, 


wobei die Koeffizienten a; Funktionen von x sind, a;=a,(x), die fiir alle ganz- 
zahligen x20 endliche Werte annehmen mégen, oder durch Umformung von 
(5.1) nach einer bekannten Formel die Gestalt 


5.2) Li=f(x-n) 4 bf (*4+n—1)+---+, f(x), 6,=6,(%),. v=1,...,%. 


Wir setzen voraus, da8 die Funktion 0,,(«) =-0 ist fiir alle ganzen x20. 


Ist dann ,(x) ein System linear unabhangiger Lésungen der homogenen 
Gleichung, Lf=0, das fiir alle ganzen x=O definiert und endlich ist, so gilt 
bekanntlich® [wegen b,, (x) ==0 («=0)] fiir die Determinante 


(53) /D(4)—det(y,(% + u—1)) =0- & w=1,...,.”) firallex=0, 
und man kann mit Hilfe der Lagrangeschen Methode der Variation der Kon- 
stanten® einen Operator L~ konstruieren, der die Gleichung L {= w (x) lost. Dieser 


Operator L“ ist sicher vollstetig im Sinne der Norm® in jedem endlichen Grund- 
intervall OX *<f (x ganz). Besitzt daher das Gleichungssystem (1.3) mit w= (x) 
fiir jeden Vektor c=(c,) ein eindeutiges Losungssystem’? /,=/(k), R=0,...,n—1, 
und ist M(f) eine stetige und beschrankte Funktion von f, so sind die Voraus- 
setzungen 1a—c des Satzes erfiillt. 


1 Dies braucht jedoch fiir den Existenzbeweis nicht durchgefiihrt zu werden. 

2 Das heiBt, solcher 7 —m Komponenten des Vektors u, bei deren Nullsetzen die 
Gleichung Lu= 0 nur die Nullésung hat. Vgl. R. ZURMUHL [13], S. 98f. 

3 Da bei Vektoren ein stetiges lineares Funktional ebenfalls ein linearer homogener 
Ausdruck der Komponenten ist, kann natiirlich hier das Gleichungssystem geschlossen 
behandelt werden. Eine solche Aufteilung wie oben ist nur vorteilhaft, wenn bei 
festen linken Seiten der ersten m Gleichungen die Losbarkeit fiir verschiedene Koet- 
fizientensysteme der iibrigen »—m Gleichungen untersucht werden soll. (1.7) ergibt 
dann einen linearen Ausdruck dieser Koeffizienten. 

4 Ak f(x) sind die absteigenden Differenzen k-ter Ordnung. Siehe etwa G. SCHULZ 
[12], S. 69 oder A. O. GELFoND [14], S. 23. 

5 Siehe etwa A. O. GELFoND [1/4], S. 266ff. 

P 

6 Zum Beispiel ||/(*)| =ZIf(A)| oder |7(-)| Seal u.a. 

7 Die stetige und eindeutige Abhangigkeit der Lésungen von den Anfangswerten 
ist hier gewdahrleistet. 
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Diese Bedingungen sind z.B. bei der Differenzengleichung 


A” f(x) = a(x) sin f(x + 2) + B(x) 


erfiillt, wenn als Fundamentalsystem etwa y,(x)=2°, P=0,...,~—1, gewahlt 
wird. 
Als Nebenbedingungen (1.6) treten Bedingungen der Gestalt auf: 


(5.4) Sur Tee 5 mit p=Ty) (%, 2 0, ganz) (u=1,..., 7) 


(opi 
mit fiir >) | f,|<0o beschrankten, stetigen Funktionen H,,. 
k=l 
Fiir die Loésbarkeit von (1.1), (5.4) [w=f(x)] ist dann die Bedingung 


N 
det («,,,) = det/ ue Cua Pv (%) = 0 
o=1 


hinreichend. Der spezielle Fall der Anfangswertaufgabe {(k) = A,,k=0,...,1—1, 
ist mit D(0)==0 (in (5.3)) noch einmal in den Bedingungen (1.7) enthalten. 


6. Bei nichtlinearen Differentialgleichungen oder Integrodifferentialgleichungen* 
bzw. Systemen solcher Gleichungen kann man, sofern ein Fundamentalsystem 
vorliegt, mit Hilfe der Methode der Variation der Konstanten? einen vollstetigen 
Operator L+ konstruieren. Es ist jedoch darauf zu achten, daB im Falle des 
Auftretens héherer Ableitungen in (1.6) die Stetigkeit der Funktionale gewahr- 
leistet bleibt. Die stetige Abhangigkeit der rechten Seite von (1.3) von den 
Konstanten c,, v=1,...,, entspricht dann der stetigen Abhangigkeit der L6- 
sungen von (1.1) von den Anfangswerten. Im Spezialfalle emes Systems von 
Differentialgleichungen 1. Ordnung 


(6.1) i BMV aes Veh eae Bee 


oder einer expliziten Differentialgleichung n-ter Ordnung 


(6.2) Ve (as Dae. weenie ae) 


stellen die Konstanten c; unmittelbar bzw. bis auf eventuelle konstante Faktoren 
die Anfangswerte y;(x)) bzw. y® (x9), p=0,...,—1, dar. 

Die Bedingungen (1.6) erfassen z.B. die schon sehr allgemeinen Randbe- 
dingungen der Form (bei einer einzelnen Gleichung (6.2) -ter Ordnung) 


N ;n-1 ; 
(63) DL Deny) tla =A,  %=0..,.0—4, 


= 


mit Funktionaloperatoren 
bik Trp T 


(6.3 a) LalyQl=f So Sei ltis Tr) ¥(%) 4%... aT 


Vk Vk Qik 


1k? 


: Im allgemeinen mu8 man dabei voraussetzen, da® die Integralanteile zu Mu 
gehoren. Operatoren L mit co-dimensionalem Nullraum erfordern eine gesonderte 
Betrachtung. 


®» Siehe etwa L. Cottarz [8], S. 83 und 96 sowie W. M. WuyBurNn [2] und die 
dort angegebene Literatur. 
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wobei die Konstanten a;,, b;,, *;, eimem Grundintervall [a, 6] angehéren!, oder 
Aufgaben mit Randbedingungen (bei einem System (6.1)) der Form 


(6.4) f dF) y() =c 
y 
mit einer Matrix F(t), deren Elemente in dem zugrunde gelegten Intervall von 
beschrankter Variation sind, wie sie R. Conti [6] behandelt. Aufgaben mit all- 
gemeineren nichtlinearen Randbedingungen der Form (1.6) sind meines Wissens 
noch nicht behandelt worden. Es ergeben sich hierfiir jedoch dieselben Bedin- 
gungen, da die Nichtlinearitat auf der rechten Seite von (1.6) auftritt. Dagegen 
sind vielfach Spezialfalle von (6.3), (6.3a) bzw. (6.4), die auch in den Anwen- 
dungen auftreten?, getrennt untersucht worden. 


Hier ergibt sich aus (1.7) z.B. mit (6.2), (6.3) mit einer in dem Streifen [a, 0], 
|y|< co beschrankten Funktion f, die so beschaffen ist, daB die Lésungen der 
Anfangswertaufgabe (%) =p, P=0, ..., m—1, stetig von den Anfangswerten 
abhangen®, unmittelbar die fiir die issues hinreichende Bedingung, da die 
Determinante 4A der Matrix 


N 


(6.5) oon = (3 > (on eet ds eal 


bh 
p! 


)): 1 Pi Oe ca HA 


mit x?,=41 fiir x;,=0 von Null verschieden ist. Dabei sind als Fundamental- 
: p-i 
system die Funktionen p, = ae p=0,...,n—1, mit den Ableitungen oY = = isa 


fiir 7 <4, yp! =0 fir; > >, gewahlt worden; und die Numerierung haben wir zweck- 
maBig von 0 bis m — 1 laufen lassen. 


Fiir den Spezialfall der Nicolettischen Aufgabe y(x,;)=A;, 1=0,...,n—1, 
erhalt man beispielsweise aus (6.3), (6.5) (mit x?=1 fiir x;=0) 


det («;») act( 3° a, act (71) eT ) £0 
et (a;,) = de se me ad Wee ok oe ee |e —. Xy— Xy 
p pa 19 (p—j)! p! tf yl in 


falls x, = x, fiir uy gilt*. 
Die Aufgabe (6.2), y' (x, =A;, 0<rSn—1 ist dagegen im allgemeinen nicht 
lésbar®: Hier ist «;)=0 fiir alle 7, also 4=0. 


1 Die Aufgabe, bei der im Argument der L;, noch Ausdriicke y'?) (x), p=1, ...,n—1, 
additiv hinzukommen, ist nicht allgemeiner, da die Ableitungen durch partielle Inte- 
gration beseitigt werden k6nnen. 

2W.M. WHYBURN [3]. 

8 Gewohnlich wird hierfiir vorausgesetzt, daB die Funktion / den Fat aes cre 
Bedingungen geniigt, d.h. fe Grandstreifen [a, b], |v? oe) =O), ==aly bbe 
festes x stetig beziiglich y), p=o, m—1, und - feste y(?) integrabel beztiglich 
w ist und ferner eine Un ten a Lipschitz-Bedingung bzw. eine schwachere 
Bedingung von McSuanéE erfiillt. Siehe etwa W. M. Wuysurn [2], R. Conti [6] und 
E. J. McSuHaneE [9], S. 349. 

4 Die bis auf die Faktoren 1/p! Vandermondsche Determinante («;») stellt in den 
x,, v=0,..., m—1, ein Polynom (m—1)-ten Grades mit den Nullstellen Ay, WY, dar. 

5 Zum Beispiel hat schon y’’=0, y’(0)=0, y/(1)=1 keine Losung. 

24* 


356 Hans EHRMANN: 


Fiir den speziellen Fall der Anfangswertaufgabe y (x9) =A;, 7=0,..., 4-1, 
wird (1.7) die Wronskische Determinante, die nach den Voraussetzungen =-0 ist, 


== a WW? ; 
und im speziellen Fall (6.2) mit dem Fundamentalsystem? pee a ergibt 


sich aus (6.5) die Einheitsmatrix «;5= jp. 


7. Von besonderem Interesse sind haufig die periodischen Lésungen y (x + X) = 
y (x) der Differentialgleichung. Hier ergibt sich aus dem Existenzsatz unmittelbar 
der Satz von R.Cont12, daB die Gleichung (1.1) [v=y(x) bzw. u= (vy; (x)) | 
unter den Voraussetzungen des Satzes mindestens eine periodische Lésung besitzt, 
wenn die Gleichung Lu—@ als einzige periodische Lésung die Nullésung hat. 
Allgemeiner kénnen wir nach einer Lésung mit einer bestimmten Eigenschaft A 
fragen. Es ergibt sich aus dem Existenzsatz dann die einfache Folgerung: 

Sind unter gewissen zusdtzlichen Voraussetzungen B fiir die Gletchung (1.1) 
gewisse Nebenbedingungen der Form (1.6) hinreichend fiir die Existenz einer Losung 
mit der Eigenschaft A, so existiert mindestens eine solche Lésung, wenn B erfiillt 
ist und die vollhomogene Aufgabe Lu=O, F,[u|=0, mit den betreffenden Bedin- 
gungen nur die Nulldsung hat. 

Der Satz von R. Conti ist jedoch in manchen Fallen nicht anwendbar zum 
Nachweis einer periodischen Lésung. Zum Beispiel gilt dies schon fiir den Fall 
(6.2) oder bei Systemen (6.1) der Gestalt y’=A(x)y+a(x, y) [A nxn-Matrix, 
y, a nm-Vektoren], wenn nach Ausmultiplikation A(x)y eme Komponente des 
Vektors y nicht auftritt. Dies ergibt sich auch aus den Bedingungen. Zum Beispiel 
erhalt man fiir das Randwertproblem 


(FAyO cyl ea (apna sassy) eel BRAS ay Syl ene vee ee!) IR A aXe 0) 


mit den fiir die Periodizitat einer L6sung notwendigen und hinreichenden Rand- 
bedingungen 


(72) y (0) = y (XJ=0, 4=0,...,%—4, aus® (6.5) 
| ofr sXe fine semtse et. 
X&ip=4 (Pt)! (p—7)! 
| 0 fir i>p 


und daher «,,—0 fiir s=n—1, p=0,...,%—1. (1.7) ist also nicht erfillt. 

Man kann jedoch auch hier die ,,Methode der Randwertaufgaben‘‘ zum Nach- 
weis periodischer Lésungen anwenden, wenn die Funktion f auBer der Periodizitat 
noch zusatzliche Eigenschaften besitzt*. Ist z.B. f so beschaffen, daB mit y= y (x) 
auch y= (x)= — p(X — x) Lésung ist, so folgt aus den Gleichungen 


(7.3) pR (a) = — g@)(X — x), p(y Spe x), S01 2... 


LY) 


* Man beachte wieder, daB wegen yy=1 in der Determinante 0°=1 zu setzen ist. 
Ke CONaI (6) tSa didi7e 
® Mit L;,=0, N=2, %1=0, %j2=X, C4j;,—=(—1)*41-6;,. 
* Dies ist bei Differentialgleichungen 2. Ordnung mehrfach durchgefiihrt worden. 
Siehe z.B. F. LETTENMEYER [19] und H. ExRMaNN [22]. In [22] wird eine voll- 
standige Ubersicht iiber diese Falle fiir Differentialgleichungen 2.Ordnung gegeben. 
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wobei die Ableitungen nach dem Argument zu nehmen sind, daB die Rand- 
bedingungen 


(7.4) ys) (0) = 0, yeo(=) =o, SSO eu 


hinreichend fiir die Periodizitat der Lésung sind}. 


Es sei nun gerade. Dann haben wir in (7.4) genau » Randbedingungen. 
Diese lassen sich erfiillen: 
Man erhalt aus (6.5) mit 
O. : fii . 
aie N=1, C =| 7,0 ue gerade 7 . 
;i-1 flr ungerade z, 
fiir gerade 1 


0) 
MN ANG ee ' (x; =1 fiir x;=0), 
=a fiir ungerade 2 


nach leichter Rechnung 


Spezielles Beispiel: 
(75) yO4 t(y)=csinwx, mngerade, X= 


25 
By ~ 

Ist f(y) eine ungerade Funktion, so ist hier die Bedingung erfiillt, daB mit 
y (x) auch w(x) = — y(X — x) Losung ist, wie man leicht durch Einsetzen erkennt. 
Ist f(y) dariiber hinaus noch stetig und beschrankt, so hat (7.5) eine periodische 
Lésung mit der Periode X. 


Zum Beispiel besitzt die verallgemeinerte Duffingsche Differentialgleichung? 
y™ + ¢, sin y = cy sin@ % 


fiir eine beliebige gerade Ordnung m mindestens eine periodische Lésung mit der 
Periode X = 27/w. 
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